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The  primary  goal  of  this  study  is  to  determine  whether  non-invasive  magnetic  resonance  (MR) 
techniques  can  distinguish  between  slow  and  rapidly  growing  and  metastatic  prostate  tumors. 
This  is  particularly  important  in  prostate  cancer  where  30%  of  men  over  the  age  of  50  have 
prostate  cancer  at  autopsy  but  only  10%  of  men  develop  prostate  cancer.  Reliable  methods  do 
not  exist  to  detennine  which  cancers  are  aggressive  and  need  to  be  treated  vs.  patients  who  could 
undergo  watchful  waiting.  A  second  goal  is  to  determine  if  anti-angiogenesis  agents  can  be  used 
as  chronic  low  toxicity  therapy  for  “newly  diagnosed”  tumors  and  as  adjuvant  to  “curative” 
therapies  and  if  MR  techniques  can  be  used  as  an  early  or  a  priori  marker  of  tumor  response. 
We  propose  that  non  invasive  measurements  of  tumor  vascular  volume,  permeability,  choline 
and  lactate  will  predict  tumor  aggressiveness  (growth  rate,  tendency  to  metastasize).  The 
hypothesis  is  based  on  data  which  indicates  that  tumor  growth  rates  and  metastases  are  related  to 
angiogenesis  which  can  be  detected  by  dynamic  contrast  enhanced  magnetic  resonance  imaging 
(DCE-MRI). 

In  addition,  we  propose  to  determine  whether  chronic  anti-angiogenic  therapy  can  be  used  in 
small  “newly  diagnosed”  tumors  and  as  an  adjuvant  to  radiation.  It  will  also  be  determined 
whether  MRSI  and  DCE-MRI  can  predict  response  to  anti-angiogenic  therapy.  We  will  also 
study  the  effect  of  anti-angiogenic  therapy  in  small  “newly  diagnosed”  tumors  and  also  as  an 
adjuvant  post  radiation  to  determine  if  it  delays  tumor  growth  and  metastases.  Tumor  doubling 
times  and  number  of  metastatic  lesions  in  the  lung  will  be  the  biological  outcome  measures. 
Developing  chronic  treatment  modalities  that  delay  (or  obviate)  the  need  for  radical  therapy  will 
enhance  treatment  of  prostate  cancer  and  possibly  quality  of  life.  Similarly,  developing  a 
treatment  that  enhances  response  to  radiation  could  potentially  enhance  the  cure  rate  or  delay 
recurrence  which  might  alter  subsequent  therapy.  The  techniques  applied  here  are  available  or 
readily  implemented  on  clinical  scanners  (choline  and  lactate  detection;  DCE-MRI)  so  the 
research  is  highly  translational. 


Body 
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Data  Presentation:  Biological 

A.  Data  acquisition  (does  not  differ  significantly  from  last  year’s  report): 

For  MR  studies,  tumor  bearing  rats  were  anesthetized  with  isoflurane  (1.0-2. 5%)  and 
compressed  air.  A  24  gauge  tail  vein  catheter  (Terumo  Medical  Corporation,  Elkton,  MD)  was 
inserted  into  the  tail  vein  prior  to  placing  the  tumor  bearing  rat  in  the  magnet.  MRS  (magnetic 
resonance  spectroscopy)  lactate  determination  and  dynamic  contrast  enhanced  magnetic 
resonance  studies  were  conducted  for  tumors  in  the  range  of  300-2860  mm3  in  the  Dunning  H 
and  200-2860  in  the  Dunning  R3327AT  model.  The  studies  were  grouped  based  on  similar 
tumor  volumes  for  the  data  analysis.  The  studies  discussed  in  the  previous  report  which  are 
briefly  resummarized  since  we  have  more  data.  For  the  Dunning  R3327-AT  tumor,  the  groups 
were  110-265  mm3  (I),  478-649  mm3  (II),  1035-1207  mm3  (III),  1409-1708  mm3  (IV)  and  1941- 
2453  mm  (V).  For  the  Dunning  H  models,  the  studies  were  similarly  grouped  as  227-259  mm 
(I),  430-546  mm3  (II),  1000-1236  mm3  (III),  1382-1715  mm3  (IV)  and  1886-2255  mm3  (V). 
Volumes  are  corrected  for  1.3  mm  skin  thickness.  Upon  completion  of  the  lactate  studies,  DCE- 
MRI  (dynamic  contrast  enhanced  magnetic  resonance  imaging)  studies  were  performed  on  the 
rats  without  removing  them  from  the  magnet.  Tumor  bearing  rats  were  typically  studied  at  2-5 
volumes.  After  the  last  study,  the  tumor  bearing  rats  were  injected  with  pimonidazole  (60mg/kg) 
and  excised  60  minutes  later.  Pimonidazole  and  H&E  staining  were  performed  on  tumor  slices  of 
8  micron  thickness.  Some  rats  were  specifically  studied  at  smaller  volumes  for  comparison  to 
the  MR  data;  in  the  latter  category,  the  study  was  terminated  at  smaller  tumor  volumes  to  obtain 
the  histology  data. 

MR  experiments:  Magnetic  resonance  experiments  (MRS)  Experiments  were  performed 
on  a  Bruker  4.7  T,  40  cm  bore  animal  scanner.  A  home-made  2  turn  volume  copper  foil  coil  with 
a  25  mm  diameter  was  used  as  a  transmit-receive  radiofrequency  coil  for  all  studies.  Scout 
images  were  obtained  in  the  three  orthogonal  directions  to  ensure  that  the  tumor  was  positioned 
within  the  field  of  view.  MRS  lactate  detection  was  performed  using  the  selective  multiple 
quantum  coherence  (Sel-MQC)  (1)  sequence,  using  selective  15  ms  single-lobe  Sine  pulses.  The 
Sel-MQC  sequence  selectively  detects  lactate  by  selecting  the  zero  quantum  (ZQ)  — >■  double 
quantum  (DQ)  coherence  transfer  pathway.  Phase  cycling  gradient  combinations  of  gl:g2:g3  = 
0:-l:2  with  duration  Si  =  S2  =  2  ms,  S3  =  4  ms,  and  amplitude  of  24  G/cm  are  used.  The  pulse 
sequence  parameters  for  the  lactate  editing  experiments  included  512  data  points,  8  averages, 
TR=2  s  and  a  spectral  width  of  2500Hz.  Two-dimensional  chemical  shift  imaging  was  performed 
with  a  matrix  size  of  16x16  FOV  of  40  mm  (2.5  x  2.5  mm  in  plane  resolution)  and  a  5  mm  slice 
thickness  generated  by  a  1ms  three-lobe  Sine  pulse.  The  2D  CSI  lactate  map  was  coregistered 
with  T2-weighted  images  of  5  mm  slice  thickness  from  the  center  of  the  tumor. 

Upon  completion  of  the  lactate  imaging,  DCE-MRI  studies  were  performed.  Based  on  the 
previous  scout  images,  3  slices  in  the  sagittal  direction  were  imaged  using  a  gradient  echo  fast 
imaging  (GEFE)  sequence  with  proton  density  weighted  parameters.  The  acquisition  parameters 
included  2  mm  slice  thickness,  a  gap  of  0.2  mm,  TR/TE/flip  angle  =  500/3/30°,  256  x  128 


matrix,  and  2  excitations  per  phase  encoding  step.  Subsequently  the  dynamic  images  were 
obtained  using  the  same  parameters  as  above  except  that  the  TR  was  reduced  to  50ms  and  288 
data  sets  were  obtained.  Gd-DTPA  (0.2  mM/kg)  was  injected  intravenously  after  the  first  5 
images  of  the  dynamic  train  had  been  acquired  and  these  images  were  used  to  measure  the 
precontrast  signal  intensity  under  saturation  conditions. 


Data  Processing:  DCE-MRI  data  were  analyzed  using  the  two  compartment  model  of  Hoffman  et 
al  (2)  to  evaluate  the  rate  constant  “kep”  and  Akep  (A  =  amplitude  -  one  of  the  parameters  from 
the  fitting  of  the  data  to  the  model;  kep  =  rate  constant  for  transfer  of  Gd-DTPA  from  the 
interstitial  compartment  to  the  plasma)  and  slope  (rate  of  increase  of  signal/time)  using  IDL 
programs  written  previously  for  clinical  studies. 


Figure  la 


Figure  1.  Tumor  volume  (mean  ±  SEM)  measurements  of 
Dunning  R3327  AT  (A)  and  Dunning  H  (B)  tumors.  Note 
different  time  scales 

then  used  to  quantify  the  tumor  voxel  lactate. 


MRSI:  Spatial  Fourier  transform 
and  superimposition  of  the  spectral  grid 
on  the  corresponding  T2-weighted 
image  were  performed  using  the  3DiCSI 
software  package  (courtesy  of  Truman 
Brown,  Ph.D.,  Columbia  University). 
The  voxels  within  the  tumor  were  then 
identified  and  the  free  induction  decay 
(FID)  from  each  tumor  voxel  was 
extracted.  The  data  for  each  FID  was 
input  to  the  jMRUI  software  package 
and  the  lactate  resonance  was  fit  in  the 
time  domain  using  AM  ARES 
(Advanced  Method  for  Accurate, 
Robust,  and  Efficient  Spectral  fitting) 
(3).  The  reference  phantom  2DCSI  set 
was  processed  similarly.  For 
quantitation  purposes,  each  tumor  voxel 
was  referenced  to  a  phantom  voxel  at 
the  exact  same  location.  The  tumor 
voxel  lactate  peak  area,  the  reference 
voxel  lactate  peak  area,  the  appropriate 
T1  saturation  factors  and  the 
comparative  coil  loading  parameters  for 
tumor  and  phantom  studies  and  were 


D.  Results 


In  last  year’s  report  we  displayed  growth  curves  for  the  Dunning  R3327AT  and  Dunning  H 
tumors,  preliminary  CSI  spectra  for  Dunning  R3327AT  and  H  tumors  (volume  range  from  about 
300  mm3  to  3000mm3),  and  dynamic  contrast  enhanced  MRI  data  on  the  same  tumors.  We 
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have  by  now  completed  the  final  analysis  and  submitted  the  data  (Appendix  2).  We  summarize 
the  results  below  (which  are  taken  from  Appendix  2). 


In-vivo  lactate  measurement 

The  anaplastic,  faster-growing  R3327-AT 
cell  line  had  an  in  vivo  tumor  doubling 
time  (TDT)  of  3.3  ±  0.17  days.  The  well- 
differentiated,  slower  growing  Dunning  H 
had  a  doubling  time  of  22.8  ±  1.95  days 
(Figure  1)  (p<0.001),  compatible  with 
previous  studies  (4). 

Figure  2 A  shows  a  series  of  T2  weighted 
MR  images  of  5mm  thickness  co¬ 
registered  with  2D-CSI  lactate  spectra 
obtained  on  the  Dunning  R3327  AT  tumor. 
Lactate  peaks  in  the  tumor  region  are 
highlighted.  Lactate  was  not  detected  in 
R3327-AT  tumors  studied  in  the  volume 

■j 

range  110-265  mm  (group  I;  data  not 
shown).  As  tumor  volumes  increased  to  the 
range  of  478-649  mm  (group  II),  the 
average  voxel  lactate  concentration 
detected  was  1.54  ±  0.54  mM  (Fig  3  A).  In 

■j 

group  III  (tumor  volume  1035-1207  mm  ), 
the  lactate  concentration  increased  to  4.79 
±  1.67  mM.  In  Groups  IV  (1409-1708 
mm3)  and  V  (1941-2453  mm3),  lactate  was 


Figure  2A.  2DCSI  lactate  spectra  of  R3327  AT  prostate 
tumor  model  as  a  function  of  tumor  volume,  1:613  mm3, 
11:1054  mm3, 111:1455  mm3  and  IV:  2207  mm3.  T2  weighted 
saggital  MR  tumor  image  is  co-registered  with  2DCS1  lactate 
spectra.  The  peak  intensity  of  the  lactate  resonance  is 
heterogenous,  even  within  a  specific  tumor.  Figure  2B. 
2DCSI  lactate  spectra  of  Dunning  H  prostate  tumor  model  as 
a  function  of  tumor  volume,  A:  246  mm3  B:512  mm3,  C:1075 
mm3,  D:1382  mm3  and  E:  2255  mm3.  T2  weighted  saggital 
MR  tumor  image  is  co-registered  with  2DCS1  lactate  spectra. 
Note  the  absence  of  lactate  in  any  of  the  voxels 


4.64  ±  1.26  mM  and  2.71  ±  0.93  mM 
respectively  (Fig  3A,  Table  1).  The 
variation  of  average  lactate  concentration 
with  growth  is  shown  graphically  in  Figure 
3A  and  demonstrates  the  initial  increase 
and  subsequent  decrease  in  lactate  noted 
with  tumor  growth.  Results  of  the  non¬ 
linear  mixed  effects  analysis  showed  very 
significant  curvature  in  the  relationship 
between  lactate  and  tumor  volume 
(p=0.008)  namely,  an  initial  increase  followed  by  a  decrease  after  reaching  a  peak  that  varied 
among  rats  at  volumes  of  1048  +/-  467  (SD)  mm  .  There  was  a  significant  quadratic  or  peak 
effect  in  the  relationship  between  lactate  per  voxel  and  tumor  volume  (p<0.001)  in  the  R3327- 
AT  tumor  with  significant  variation  among  rats  in  the  volume  at  which  the  maximum  expected 
value  of  (log)  lactate  occurred. 


Figure  3A.  Measurement  of  average 
lactate  in  R3327  AT  tumor  model  as 
a  function  of  tumor  volume.  Tumors 
in  the  volume  range  of  1 10-  265 
mm3  did  not  show  lactate.  Figure 
3B.  Evaluation  of  tumor  core/rim 
lactate  of  R3327  AT  tumor  model  as 
function  of  tumor  model  (solid 
fdled  bar  plots=  core  lactate;  empty 
bar  plots=  rim  lactate).  Lactate 
initially  increases,  shows  a  plateau, 
and  then  decreases  in  the  whole 
tumor,  rim,  and  core 
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When  the  lactate  concentration  in  the  outer  voxels  of  the 
tumor  (“rim”)  was  compared  to  all  the  other  voxels  (“core”) 
(Fig.  3B),  lactate  was  significantly  higher  in  the  core 
compared  to  the  rim,  independent  of  tumor  volume  (p<0.001). 
The  lactate  concentration  vs  tumor  volume  curve  for  the  core, 
rim  and  whole  tumor  increases  to  a  maximum  followed  by  a 
decrease  after  reaching  a  peak  or  plateau  in  groups  III  and  IV. 
This  was  characterized  by  a  quadratic  term  in  the  statistical 
model,  with  a  coefficent  that  was  significantly  less  than  zero, 
reflecting  a  downturn  in  lactate  concentration  with  volume 
after  a  peak.  Inclusion  of  an  additional  parameter  in  the 
statistical  model  indicated  weak  evidence  (p=0.09)  that  on 
average,  the  lactate  maximum  concentration  occurred  at 
larger  tumor  volumes  for  the  core  voxels  compared  to  the  rim. 


A  representative  set  of  2DCSI  spectra  of  Dunning  H  prostate 
model  at  different  volumes  is  represented  in  Figure  2B 
(groups  I  to  V).  Lactate  was  not  detected  in  any  tumor, 
regardless  of  the  volume. 

DCE-MRI  analysis 
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Dunning  R3327-AT 

The  average  Akep  for  group  I  tumors  was  1.59  ±  0.53  /min. 
With  an  increase  in  tumor  size  (478  -  649  mm  ),  Akep 
decreased  to  0.65  ±  0. 12.  With  further  increases  in  tumor 
size,  the  Akep  value  for  group  III  and  IV  tumors  were  0.55  ± 
0.1 1  and  0.36  ±  0.04  respectively.  The  largest  tumors  (group 
V)  had  a  value  of  Akep  of  0.34  ±  0.05  (Fig.  4A;  Table  2A). 
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We  analyzed  whether  the  rim  and  core  Akep  measurements 
differed.  The  same  criterion  for  classifying  the  rim  and  core 
voxel  was  adopted  as  was  used  in  the  lactate  analysis  (Fig. 
4B).  The  average  tumor  rim  and  core  Akep  value  of  the 


lowest  tumor  volume  cohort  were  1.61  ±  0.52  /min  vs.  and 
0.15  ±  0.14  /min  respectively,  but  only  one  of  the  4  animals 
had  a  tumor  large  enough  to  delineate  core  voxels  so  the 
biological  significance  of  this  difference  is  unclear.  With 
the  increase  in  tumor  volume  to  478-649  mm  ,  the  tumor  rim 
and  core  Akep  values  were  0.71  ±  0.12  /min  vs.  0.27  ±  0.07 
/min.  As  expected,  there  are  a  greater  number  of  voxels  in 
the  rim  than  in  the  core.  At  higher  tumor  volumes  (1035-1207  mm  and  1409-1708  mm  ),  the 


Figure  4 A.  Whole  tumor  average  Akep 
(/min)  of  Dunning  R3327  AT  tumor  model 
as  a  function  of  tumor  volume.  Akep 
decreases  with  tumor  volume.  Fig  4B. 
Evaluation  of  tumor  core/rim  Akep  (/min) 
of  Dunning  R3327  AT  tumor  model  at 
different  tumor  volumes. 


ritn  and  core  Akep  values  decreased  to  0.63  ±  0.12  /min  vs  0.36  ±  0.12  /min  and  0.50  ±  0.07  /min 
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vs  0.15  ±  0.03  /min  respectively.  At  the  largest  tumor  volumes  the  rim  and  core  Akep  values 
were  0.47  ±  .07/min  and  0.14  ±  0.02/min  respectively.  The  rim  Akep  values  were  significantly 
higher  than  core  values  for  the  Dunning  R3327-AT  tumors  (p<0.001) 

DH  tumors 

Slow  growing  Group  I  Dunning  H  tumor  had  an  average 
Akep  value  of  1.23  ±  0.12/min  (Fig.  5A,  Table  2B).  This 
decreased  to  1.06  ±  0.23/min,  0.84  ±  0.21/min  and  0.86  ± 
0.08/min  in  groups  II,  III  and  IV  tumors  respectively.  With 
continued  tumor  growth,  the  whole  tumor  Akep  values 
decreased  (Group  V  =  0.46  ±  0.04/min).  We  subsequently 
subdivided  the  analysis  to  study  the  rim  and  core  as  we  had 
done  previously  for  the  Dunning  R3327-AT  tumor  except 
in  the  Group  I  tumors  wherein  there  were  no  core  voxels 
because  of  the  small  tumor  size.  The  rim  and  core  Akep 
values  both  decreased  with  growth  (p<0.001). 

The  Dunning  H  exhibited  heterogenous  perfusion  as 
manifested  by  significant  differences  in  Akep  value 
between  the  rim  and  core  (p<0.023).  The  Dunning  H 
tumors  had  significantly  higher  Akep  values  than  the 
R3327-AT  tumors,  as  measured  in  the  whole  tumor 
(p<0.001),  the  core  (p<0.001)  and  the  rim  (p<0.001).  Akep 
declined  progressively  over  the  range  of  tumor  sizes  in  the  Dunning  H  tumors  (p<0.001), 
particularly  at  the  largest  tumor  volumes.  After  applying  an  appropriate  normalizing 
transformation  to  Akep  (e.g.,  cube  root),  the  rate  of  decrease  was  linear.  This  decrease  was  also 
noted  if  the  data  were  analyzed  separately  for  the  rim  and  core  (p<0.001). 

Histology  and  Immunohistochemistry  Correlation 

The  Dunning  R3327-AT  tumors  have  necrosis  present  even  at  small  tumor  volumes  (Fig.  6a,  b, 
Fig.  7A).  In  contrast,  the  Dunning  H  tumors  had  no  observable  necrosis  even  at  volumes 
approaching  2300  mm  ,  (Figure  6c).  In  the  R33327-AT  tumor,  we  found  that  with  increasing 
tumor  volume,  the  necrotic  fraction  increased  significantly  (R  =  0.808  ;  p  <0.001)  (Fig,  7A).  A 
linear  correlation  between  the  concentration  of  lactate  measured  and  the  viable  tumor  volume 
was  noted  indicating  that  the  concentration  of  lactate  increases  with  tumor  growth  corrected  for 
viability  (R2  =  0.764;  p=0.005).  (Figure  7B) 

Tumor  sections  sequential  to  those  used  for  histological  staining  were  assessed  for  the 
distribution  of  the  hypoxia  tracer  pimonidazole  by  immunofluorescence  staining  (Fig.  6(d),  (e), 
(f)).  We  observed  that  the  fractional  pimonidazole-positive  tumor  area  increased  with  tumor 
volume  (R  =  0.720,  p  =  0.007).  Comparison  of  pimonidazole-positive  tumor  fraction  to  lactate 
values  also  showed  a  statistically  significant  positive  correlation  (R  =  0.602,  p  =  0.023;  Fig. 
1C). 
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Figure  5 A.  Whole  tumor  average  Akep 
(/min)  of  Dunning  H  tumor  as  a 
function  of  tumor  volume.  Figure  5B. 
Evaluation  of  tumor  core/rim  Akep 
(/min)  of  Dunning  H  tumor  model  over 
the  volume  range  studied. 
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SS-Sel-MQC  (Spectral  Selective  MQC) 

In  the  previous  report,  we  discussed  a  novel  modification  of  the  SelMQC  sequence  using 
binomial  spectral-selective  pulses  (SS-SelMQC).  Frequency  selective  excitation  pulses  were 
employed  with  suitable  phase  cycling  of  a  binomial  sequence  -a-u/4)^  to  selectively  excite 

either  lactate  methyl  CH3  resonances  c|)2=x,  -x)  or  methylene  CH  resonances  (cj)i,  (|>2=x,  x). 
Frequency  selective  inversion  was  achieved  using  [(^2)t-a-(^/2)_J  for  the  lactate  methyl  CH3 
resonances.  Chemical  shift  selection  was  achieved  by  adjusting  the  interpulse  delay  equal  to  the 
inverse  of  twice  the  difference  in  the  center  frequencies  of  maximum  and  null  excitation  bands 
(Fig.  9).  In  a  lipid  enriched  environment,  the  modified  pulse  sequence  yielded  enhanced  lactate 
signal  of  200-300%  compared  to  SelMQC.  Non-localized  proton  spectra  and  2D  CSI-lactate 
images  were  obtained  from  30mM  lactate/lipid  phantoms  and  in-vivo  R3327  prostate  animal 
tumors.  We  concluded  that  lactate  signals  detected  by  SS-SelMQC  have  a  2-3  times  higher 
signal  to  noise  compared  to  conventional  SelMQC.  The  sequence  is  currently  being  implemented 
for  lactate  detection  from  a  volume  of  interest,  which  will  be  used  to  image  the  internal  organs  of 
transgenic  animal  tumor  models.  We  have  explored  this  sequence  further  in  the  recent  year  to  try 
to  explain  the  increased  signal  to  noise. 


Figure  7A.  A  graph  of  percentage  of  necrosis  as  a 
function  of  tumor  volume  demonstrates  a  significant 
linear  correlation  between  these  two  parameters. 

Figure  7B.  A  graph  of  the  variation  of  tumor  lactate  vs. 
viable  tumor  volume  demonstrating  a  significant  linear 
correlation.  Lactate  is  noted  to  increase  with  tumor 
growth,  when  corrected  by  eliminating  the  necrotic 
tissue 

Figure  1C.  Graph  of  pimonidazole  positive  fraction  as  a 
function  of  measured  lactate  showing  a  linear 
correlation  between  pimonidazole  and  lactate 
concentration.  The  data  suggest  that  lactate  is  a 
potential  surrogate  for  measuring  tumor  hypoxia  in  this 
model. 
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Table  1  Average  lactate  concentration  and  numbers  of  CSI  voxels  in  the  Dunning  R-3327-AT  tumor _ 

Whole  tumor  Core  Rim 

Tumor  Tumor  volume  Lactate  Avg  #  Lactate  Avg  #  Lactate  Avg  # 

Group  (mm3)  (mM/L)  voxels  (mM/L)  voxels  (mM/L)  voxels 

I  110-265  0  7  0  1  0  6 

II  478-649  1.53  13  1.79  3  1.46  10 

III  1035-1207  4.78  17  6.95  6  3.68  11 

IV  1409-1708  4.64  20  7.34  8  3.02  12 

V  1941-2453  2.71  24  4.78  10  1.40  14 
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Average  Akep  values  of  Dunning  R  3327  AT  tumor 

Tumor 

volume 

(mm3) 

Tumor 

group 

Whole 

tumor 

Akep 

Rim 

Akep 

Average 
No.  of 
rim  CSI 
voxels 

Core 

Akep 

Average 
no.  of 
core  CSI 
voxels 

110-265 

i 

1.59±0.53 

1.61±0.52 

6 

0.15±0.14 

1 

478-649 

ii 

0.65±0.12 

0.71±0.12 

10 

0.27±0.07 

3 

1035-1207 

m 

0.55±0.11 

0.63±0.12 

11 

0.36±0.12 

6 

1409-1708 

IV 

0.36±0.04 

0.50±0.07 

12 

0.15±0.03 

8 

1941-2453 

V 

0.34±0.05 

0.47±0.07 

14 

0.14±0.02 

10 

Table  2B 


Average  Akep  values  of  Dunning  H  tumors 

Tumor 

volume 

(mm3) 

Tumor 

group 

Whole 

tumor 

Akep 

Rim 

Akep 

Average 
no.  of 
Rim  CSI 
voxels 

Core 

Akep 

Average 
no.  of 
Core  CSI 
voxels 

227-259 

i 

1.23±0.12 

1.23±0.12 

7 

— 

0 

430-546 

ii 

1.06±0.23 

1.07±0.23 

9 

0.76±0.27 

1 

1000-1236 

in 

0.84±0.21 

0.85±0.21 

11 

0.77±0.23 

3 

1382-1715 

IV 

0.86±.08 

0.90±0.04 

13 

0.75±0.19 

7 

1886-2255 

V 

0.46±.04 

0.53±0.05 

12 

0.38±.06 

11 

SS-Sel-MQC  (Spectral  Selective  MQC) 
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In  the  previous  report,  we  discussed  a  novel  modification  of  the  SelMQC  sequence  using 
binomial  spectral-selective  pulses  (SS-SelMQC)  (5).  Frequency  selective  excitation  pulses  were 
employed  with  suitable  phase  cycling  of  a  binomial  sequence  to  selectively  excite 

either  lactate  methyl  CH3  resonances  ((fq,  c|)2=x,  -x)  or  methylene  CH  resonances  (cf>i,  (|>2=x,  x). 
Frequency  selective  inversion  was  achieved  using  |^2)x-a-(^2)_J  for  the  lactate  methyl  CH3 
resonances.  Chemical  shift  selection  was  achieved  by  adjusting  the  interpulse  delay  equal  to  the 
inverse  of  twice  the  difference  in  the  center  frequencies  of  maximum  and  null  excitation  bands 
(Fig.  9).  In  a  lipid  enriched  environment,  the  modified  pulse  sequence  yielded  enhanced  lactate 
signal  of  200-300%  compared  to  SelMQC.  Non-localized  proton  spectra  and  2D  CSI-lactate 
images  were  obtained  from  30mM  lactate/lipid  phantoms  and  in-vivo  R3327  prostate  animal 
tumors.  We  concluded  that  lactate  signals  detected  by  SS-SelMQC  have  a  2-3  times  higher 
signal  to  noise  compared  to  conventional  SelMQC.  The  sequence  is  currently  being  implemented 
for  lactate  detection  from  a  volume  of  interest,  which  will  be  used  to  image  the  internal  organs  of 
transgenic  animal  tumor  models.  We  have  explored  this  sequence  further  in  the  recent  year  to  try 
to  explain  the  increased  signal  to  noise. 

Reasons  for  Signal  enhancement  in  SS-SelMQC 

Reduced  effects  of  J-scalar  coupling  evolution,  molecular  diffusion  as  well  as  77  relaxation  loss 
in  SS-SelMQC  are  the  likely  causes  of  the  signal  enhancement  observed  experimentally, 
compared  to  SelMQC.  These  effects  were  simulated  using  computer  simulations  of  un-localized 
lactate  signals  assuming  an  IS  spin  system.  Lac  editing  efficiency  in  SelMQC  sequence  depends 
on  frequency  selective  excitation  or  inversion  profiles  which  necessitates  the  use  of  long  pulses 
which  increased  the  MQ  evolution  period  and  affects  the  lactate  editing  efficiency  (1).  In  SS- 
SelMQC,  we  replaced  all  the  long  frequency-selective  pulses  with  short  binomial  composite  RF 
pulses  for  frequency  selection  resulted  in  shorter  MQ  evolution  delays.  This  effect  increased  the 
Lac  editing  efficiency  by  reducing  J-coupling  evolution  as  well  as  molecular  diffusion  effects  as 
described  below. 

Scalar  coupling  effects 

Approximately  24  %  of  the  increase  in  the  lactate  signal  observed  in  SS-SelMQC  (Fig.8A),  in 
comparison  with  SelMQC  (Fig.  8B)  (assuming  no  RF  pulse  imperfections),  can  be  explained 
from  increased  scalar  coupling  evolution  during  RF  pulse  finite  widths.  We  denoted  RF  pulse 
imperfection  ( W//W /J),  where  Wj  and  Wj°  are  actual  and  ideal  RF  field  strength  experienced  by 
samples  within  RF  coils.  Hence  (W/Wj0)  =  1  corresponds  to  no  RF  pulse  imperfections  and 
(Wj/Wj°)  =  0.9,  0.8  corresponds  to  the  presence  of  10%  and  20%  RF  pulse  imperfections 
respectively.  In  Fig.8A  and  Fig.8B,  the  x-axis  is  the  chemical  shift  scale  where  the  lactate 
doublet  centered  on  the  experimental  frequency  of  CH3  was  plotted  without  (  J^0)  and  with 
(J=0)  scalar  coupling  evolution  considered  during  RF  finite  pulse  widths  in  the  presence  of  RF 
pulse  imperfections.  The  y-axis  denotes  the  Fourier  transfonn  of  total  transverse  magnetization 
‘Mxy’  following  the  SS-SelMQC  and  SelMQC  sequences.  Lac  signal  intensity  was  normalized  as 
100%  when  there  were  no  scalar  coupling  effects  and  RF  imperfections.  Introducing  scalar 
evolution  during  finite  pulse  widths,  Lac  signal  intensity  drops  to  99%  in  SS-SelMQC,  further 
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FIG.  8.  Computer  simulated  lactate  spectra  generated  as  a  function  of  RF 
pulse  imperfections  and  scalar  coupling  evolutions  (J  _  0)  during  RF 
pulses  using  (a)  SS-SelMQC  and  (b)  SelMQC  sequence.  Parameters 
used  for  generating  these  line  plots  are  marked  on  the  figures.  Y-axis 
represents  percentage  lactate  signal  normalized  to  100%  when  there 
were  no  scalar  coupling  effects  (J  _  0)  or  RF  imperfections.  These  graphs 
show  simulated  spectra  of  lactate  CH3  resonance  corresponding  to  the 
4.7T  experimental  conditions  and  assuming  IS  spin  system  for  lactate 

/PH  =1  and  PH 3=  5^  with  q nailer  r.nunlinn  .//.9  7  Hy 


reduces  to  87%  with  10% 
and  59%  with  20%  RF  pulse 
imperfections.  Similarly,  in 
SelMQC,  Lac  signal 
intensity  drops  to  75%, 
further  reduces  to  62%  with 
10%  and  40%  with  20%  RF 
pulse  imperfections.  We  can 
further  minimize  signal  loss 
observed  in  simulations  due 
to  RF  imperfections  by 
customized  RF  coils 
generating  homogeneous  RF 
fields  does  this  belong  here. 
These  simulations  yielded  a 
21  %  signal  gain  for  SS- 
SelMQC  compared  with  SS- 
SelMQC.  It  should  be  noted 
that  for  patient  studies  at  the 
clinical  magnetic  field 
strengths  of  1.5T  or  3.0T, 
the  signal  loss  due  to  scalar 
coupling  effects  will  be 
more  severe  using  SelMQC 
compared  to  SS-SelMQC. 
The  reason  for  this  is  that 
frequency  selective  pulses 
used  in  SelMQC  need  to  be 
very  long  as  chemical  shift 
dispersion  is  very  low  and 
results  in  longer  tl  MQ 


evolution  delays  leading  to  higher  signal  loss.  Relatively  much  shorter  Tl  values  can  be 
maintained  by  binomial  selective  pulses  in  SS-SelMQC  sequence  resulting  in  minimum  scalar 
coupling  losses  increases  Lac  signal  sensitivity. 


T2  Relaxation  Losses 

Relatively  short  73  values  of  lactate  have  been  reported  in  in-vivo,  for  eg.,  rat  C6-glioma  with 
200ms  (6)  and  mice  tumors  (7)  with  MCa  (68  ms)  and  Colon-38  tumors  (117ms),  which  would 
lead  to  lower  Lac  detection  efficiency.  Since  the  in  vivo  T2  relaxation  times  are  short,  an 
optimal  sequence  minimizing  T2  losses  is  critical  for  accurate  measurements  and  for  absolute 
quantification  of  tissue  lactate  concentrations.  This  effect  may  be  particularly  important  in 
studies  with  small  tumor  volumes,  where  lactate  concentrations  may  be  low,  or  if  used  for  early 
cancer  diagnosis.  Since  the  total  duration  of  RF  pulses  in  the  modified  sequence  are  shorter  and 
approximately  45  ms  less  than  SelMQC  sequence,  this  enhances  the  lactate  signal  by  20% 
(assuming  T2=200ms)  in  the  SS-SelMQC  method  by  further  reducing  T2  relaxation  losses. 
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Molecular  Diffusion  Effects 

When  two  gradients  were  applied  during  MQ  evolution  period  to  choose  ZQ^  DQ 
pathways,  Lac  signal  attenuation  occurs  due  to  molecular  diffusion  of  molecules  during  the 
period  (28)  between  these  gradients  of  duration  8  ms,  following  the  Stejskal-Tanner  equation 
(35).  We  assumed  phantom  conditions  in  our  calculations  using  the  diffusion  coefficient  D=  1  x 
1CT  cm  s'  at  room  temperature.  Using  the  equation  reported  earlier  (1),  the  signal  loss  due  to 
molecular  diffusion  is  as  follows, 

5  /  S(  0)  =  cxp[-D{p(2  /  nySg)} 2  (A  -  i-  8) 

Where A  ' 1  '  A  +  <Ie ,  p  is  the  coherence  order,  de  is  gradient  dead  time,  ^  is  gyromagnetic 
ratio=2.675e8  s'1!'1’  for  our  present  SelMQC  parameters,  the  lactate  signal  loss  will  be  30  %  and 
signal  loss  in  SS-SelMQC  will  be  only  6  %.  Hence  a  signal  gain  of  24  %  was  obtained  in  SS- 
SelMQC.  The  molecular  diffusion  losses  will  be  severe  in  choosing  DQ->ZQ  pathway  to  obtain 
full  signal  intensity  of  lactate  or  at  low  magnetic  field  strengths  where  long  pulses  result  in 
higher  signal  loss. 

Anti-Angiogenesis  Therapy  and  Predicting  Response  (Works  in  Progress). 

We  are  in  the  process  of  performing  this  study  where  we  are  evaluating  a  novel  anti- 
angiogenic  drug  (R00281501)  and  this  represents  Work  in  Progress.  The  tumor  bearing  rats 

"3  "> 

were  divided  into  two  cohorts  with  similar  volumes  (trcatcd=  940  mm  vs.  control  =1031  mm 

respectively).  After  24hr  of  treatment,  there 
was  an  18  %  tumor  growth  inhibition 
(compared  to  control)  and  40  %  by  day  7 
(compared  to  the  control  group).  Control 
tumors  continued  to  grow  beyond  day  7  and 
became  necrotic  before  day  14,  hence  these 
tumor  bearing  animals  were  sacrificed.  The 
tumor  doubling  time  (TDT)  for  the  control 
group  was  3.3  +/-0.17  days,  whereas  the 
TDT  increased  to  approximately  7-8  days  in 
the  treated  cohort. 


Two  dimensional  chemical  shift  imaging  (2DCSI)  lactate  spectra  were  also  obtained  at  prior 
to  treatment  (day  0)  and  post  (24  hr,  day  7  and  day  14)  treatment.  Lactate  was  detected  prior  to 
treatment  (day  0;  5.6  +/-  1.13  mM/L(mean  +/-SEM))  and  is  significantly  (p  <  0.005)  decreased 
at  24  hrs  of  treatment  (Fig.  9  (2.9  +/-  0.78).  At  24  hr  post  treatment,  lactate  has  poor  signal  to 
noise  (SNR)  (Fig.  10).  By  day  7,  there  was  no  difference  between  the  lactate  values  measured 
compared  to  day  0.  At  the  highest  tumor  volume  (day  14),  the  tumor  lactate  decreased  to  2.6  +/- 
0.61  (p  <  0.03)  compared  to  day  0.  This  represents  preliminary  data  and  is  ongoing. 
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KEY  RESEARCH  ACCOMPLISHMENTS: 


1 .  Differences  between  slow  and  fast  growing  tumors  in  lactate  concentration  measured. 
This  has  been  validated  by  comparison  with  histology  and  pimonidazole  measurements 
suggesting  changes  in  perfusion  are  the  cause  of  lactate  detection  in  the  faster  growing 
tumors  as  they  get  larger  and  also  induce  hypoxia. 

2.  Development  of  new,  quantitative  and  more  sensitive  pulse  sequence  for  lactate 
measurements  implemented. 

3.  Use  of  lactate  measurements  as  an  early  predictor  of  tumor  response  to  anti-angiogenic 
therapy  -  ongoing 

4.  Use  of  dynamic  contrast  enhanced  MRI  to  predict  response  to  anti-angiogenic  therapy  - 
ongoing 
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REPORTABLE  OUTCOMES: 

Two  abstracts  were  selected  for  talks  in  the  previous  International  Society  of  Magnetic 
Resonance  in  Medicine  meeting.  One  of  these  has  now  been  accepted  as  a  complete  manuscript 
for  Magnetic  Resonance  in  Medicine  and  is  In  Press  (Thakur  et  al).  The  second  (Yaligar  et  al) 
has  been  written  and  submitted.  The  third  is  in  the  process  of  completing  the  data  analysis  and 
will  be  complete  in  about  30  days. 

Manuscripts 

Thakur  SB,  Yaligar  J,  and  Koutcher  JA.  In-Vivo  Lactate  Signal  Enhancement  Using  Binomial 
Spectral-Selective  pulses  in  SELective  MQ  Coherence  (SS-SelMQC)  Spectroscopy.  Magnetic 
Resonance  in  Medicine.  In  Press. 

Yaligar  J,  Thakur  SB,  Carlin,  S,  Thaler  HT,  Lupu,  ME,  Wang  Y,  Matei  CC,  Zakian  KL  and 
Koutcher  JA  Lactate  Imaging  and  DCE-MRI  as  surrogate  markers  of  Prostate  Tumor 
Aggressiveness.  Submitted 


Abstracts 

1 .  Comparative  study  of  tumor  lactate  and  tumor  vasculature  in  aggressive  and  indolent 
prostate  cancer  animal  models  by  2D-MR  Spectroscopic  Imaging  and  DCE-MRI  -  Yaligar  J, 
Thakur  SB,  Lupu,  ME,  Wang  Y,  Matei  CC,  Zakian  KL  and  Koutcher  JA.  Int  Society  of 
Magnetic  Resonance  in  Medicine  -  16th  Annual  Meeting,  Toronto,  May  2008,  Page  369 

2.  In-Vivo  Lactate  Detection  Using  Selective  MQ  Coherence  Spectroscpy:  Signal 
Enhancement  Using  Spectral-Selective  Binomial  RF  pulses  (SS-SelMQC)  -  Thakur  SB,  Yaligar 
J  and  Koutcher  JA.  Int  Society  of  Magnetic  Resonance  in  Medicine  -  16th  Annual  Meeting, 
Toronto,  May  2008, 

3.  Preclinical  evaluation  of  anti-angiogenic  agent  RO281501  on  R3327  AT  prostate  model 
lactate  MRS  and  DCE-MRI.  Yaligar  J,  Thakur  SB,  Coman  M,  Lupu  ME,  Wang  Y,  Kolinsky  K, 
Higgins  B,  Zakian  KL,  and  Koutcher  JA.  Int  Society  of  Magnetic  Resonance  in  Medicine  -  16th 
Annual  Meeting,  Toronto,  May  2008, 
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CONCLUSION:  We  have  shown  that  lactate  is  present  in  variable  amounts  in  a  very 
aggressive  tumor  model  (R3327AT)  which  is  known  to  be  radiation  resistant.  At  small  tumor 
volumes,  where  perfusion  was  better,  lactate  was  not  detected.  Conversely,  in  a  better  perfused, 
radiation  sensitive  tumor  model  (Dunning  H),  lactate  cannot  be  detected.  These  data  suggest  that 
lactate  might  be  a  marker  for  radiation  sensitivity. 

More  recently,  we  have  shown  that  changes  in  lactate  occur  in  response  to  successful  treatment 
with  an  anti-angiogenic  therapy  which  is  under  development.  This  occurs  within  24  hours  of 
treatment.  Further  analysis  of  the  data  stratifying  for  location  is  ongoing. 

We  have  developed  a  modification  to  the  pulse  sequence  for  detecting  lactate  which  is  more 
sensitive  and  hope  to  apply  it  in  future  experiments.  Our  preliminary  studies  suggest  that  it  is 
approximately  twice  as  sensitive  as  Sel-MQC  and  we  are  currently  evaluating  the  causes  for  this 
enhanced  sensitivity. 

Future  Plans 

A  major  hindrance  to  this  study  was  that  the  slow  growing  tumors  (Dunning  H  and  HI)  were 
difficult  to  grow.  Currently  the  Dunning  HI  is  growing  well  and  is  ready  for  transplanting  to  a 
larger  cohort  of  rats.  Additionally,  the  Dunning  H  was  obtained  from  a  second  source  (Dr.  Ralph 
Mason)  and  is  growing,  albeit  slowly  as  expected.  We  will  evaluate  these  lines  for  1)  presence 
of  lactate  in  the  HI,  2)  response  to  anti-angiogenic  therapy  and  the  utility  of  DCE-MRI  and 
lactate  changes  to  predict  response,  and  compare  the  lactate  measurements  with  a  PET  tracer  of 
hypoxia  (18F-fluoromisonidazole)  and  metabolism  (18FDG)to  detennine  if  lactate  could  be 
considered  as  a  surrogate  of  hypoxia  and  if  it  provides  equivalent  data  to  FDG  since  both  are 
probes  of  tumor  glycolysis.  By  the  use  of  the  new  pulse  sequence  developed  by  Dr.  Thakur,  we 
intend  to  improve  the  resolution  of  the  lactate  study  will  be  improved  to  that  of  the  PET  study 
(-2  x  2  x  2  mm).  We  have  had  extensive  experience  doing  multi-modality  studies  using 
different  instruments  and  registering  the  data  and  do  not  anticipate  any  problems  with  image 
registration  (8-10). 
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In  Vivo  Lactate  Signal  Enhancement  Using  Binomial 
Spectral-Selective  Pulses  in  Selective  MQ  Coherence  (SS- 
SelMQC)  Spectroscopy 

S.B.  Thakur,1’2*  J.  Yaligar,1  and  J.A.  Koutcher1-3 


Tumor  vasculature  and  tissue  oxygen  pressure  can  influence 
tumor  growth,  metastases,  and  patient  survival.  Elevated  levels 
of  lactate  may  be  observed  during  the  process  of  aggressive 
tumor  development  accompanied  by  angiogenesis  (the  evolu¬ 
tion  of  the  microenvironment).  The  noninvasive  MR  detection  of 
lactate  in  tumor  tissues  as  a  potential  biomarker  is  difficult  due 
to  the  presence  of  co-resonating  lipids  that  are  present  at  high 
concentrations.  Methods  were  previously  reported  for  lactate 
editing  using  the  SELective  Multiple  Quantum  Coherence 
(SelMQC)  method.  Here  we  report  a  sequence  “SS-SelMQC,” 
Spectral-Selective  SelMQC,  which  is  a  modified  version  of 
SelMQC  using  binomial  pulses.  Binomial  pulses  were  employed 
in  this  editing  sequence  for  frequency  excitation  or  inversion  of 
selective  lactate  resonances.  Lactate  detection  has  been  dem¬ 
onstrated  using  SS-SelMQC,  both  in  vitro  (30  mM  iactate/H20 
doped  with  25  pM  Gd-DTPA)  and  in  vivo  (Dunning  R3337-AT 
prostate  tumors),  and  compared  to  similar  measurements 
made  with  SelMQC.  Lactate  areas  were  measured  from  nonlo- 
calized  spectra,  one-dimensional  (ID)  localized  spectra,  and 
two-dimensional  chemical  shift  images  (CSI)  of  the  localized 
slice.  In  data  from  whole  phantoms,  the  modified  pulse  se¬ 
quence  yielded  enhancement  of  the  lactate  signal  of  2.4  ±  0.40 
times  compared  to  SelMQC.  Similar  in  vivo  lactate  signal  en¬ 
hancement  of  2.3  ±  0.24  times  was  observed  in  ID  slice-local¬ 
ized  experiment.  Magn  Reson  Med  62:000-000,  2009.  ©  2009 
Wiley-Liss,  Inc. 

Key  words:  lactate  editing;  single-scan  lipid  suppression;  spec¬ 
tral-selective  binomial  pulses;  2D  chemical  shift  imaging 

In  vivo  proton  MR  spectroscopy  (MRS)  provides  a  nonin¬ 
vasive,  biochemical  measure  of  metabolism  which  can  be 
used  to  differentiate  healthy  tissues  from  tumor  tissues 
based  on  the  metabolic  abnormalities  (1-3).  Metabolites 
may  be  a  priori  markers  of  prognosis  or  metabolic  changes 
can  be  very  sensitive  early  indicators  of  treatment  re¬ 
sponse.  As  such,  metabolic  approaches  show  promise  for 
monitoring,  and  perhaps  even  ultimately  selecting  treat¬ 
ments.  As  tumors  are  heterogeneous,  multivoxel  MRS  im¬ 
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aging  (MRSI)  provides  a  method  to  investigate  therapy 
response  and  optimization  of  individual  treatment  regime 
based  on  various  metabolite  levels  or  concentrations 
within  each  voxel,  which  might  signify  aggressive  (or  in¬ 
dolent)  disease  (4-8)  in  local  regions  of  the  tumor.  Among 
the  many  detectable  metabolites,  lactate  (Lac)  is  an  impor¬ 
tant  compound  that  reflects  elevated  tumor  glycolysis 
and/or  poor  tissue  perfusion,  which  may  result  in  acceler¬ 
ated  tumor  growth,  malignant  transformation,  and  metas¬ 
tases  (9).  Lac  is  the  endproduct  of  glycolysis  and  the 
methyl  resonance  peak  is  typically  a  doublet  situated  at 
1.3  ppm.  Different  levels  of  tumor  lactate  in  the  MR  spec¬ 
trum  were  reported  in  the  literature  (10-15).  High  Lac  was 
observed  in  extracts  from  biopsy  specimens  of  breast  tu¬ 
mors  (10,16)  and  prostate  tumors  (11).  Low  extracellular 
pH  and  high  Lac  levels  were  shown  to  be  indicators  of 
metastatic  risk  in  breast  cancer  xenografts  (12,13).  Ele¬ 
vated  Lac  content  in  biopsy  samples  was  shown  to  corre¬ 
late  with  increased  risk  of  metastasis  and  poor  patient 
survival  in  head  and  neck  cancer  (14)  and  cervical  cancer 
(15),  while  a  decrease  in  steady-state  tumor  Lac  levels 
related  to  tumor  response  to  radiation  (17,18)  and  chemo¬ 
therapy  (19,20).  It  has  also  been  shown  by  Pavel  et  al.  (21)  AQ:  i 
that  Lac  is  likely  to  be  found  in  malignant  breast  tumors 
that  may  be  a  marker  for  tumor  diagnosis.  Therefore,  non- 
invasively  measured  Lac  may  be  an  additional  character¬ 
istic  metabolic  marker  for  cancer  studies  and  may  improve 
diagnostic  specificity,  serve  as  an  early  marker  of  tumor 
response,  and  provide  functional  information  about  prog¬ 
nosis. 

Reliable  measurement  and  interpretation  of  MRS-de- 
rived  metabolites  requires  knowledge  of  the  factors  affect¬ 
ing  Lac  resonance  detection  and  quantitation.  Detection  of 
relatively  low  Lac  levels  in  tumor  tissues  using  conven¬ 
tional  MRS  techniques  such  as  PRESS  (22)  and  STEAM 

(23)  is  a  challenging  task  due  to  the  presence  of  the  intense 
co-resonant  lipid  (Lip)  signals.  To  differentiate  Lac  from 
the  co-resonating  Lip  resonances,  spectral  editing  tech¬ 
niques,  such  as  spin-echo  (SE)  J-difference  spectroscopy 

(24)  and  multiple-quantum  (MQ)  coherence  filters  (25-28), 
have  frequently  been  employed.  Although  the  J-difference 
approach  detects  lactate  with  100%  signal  sensitivity,  the 
subtraction  of  sequential  scans  with  different  echo  times 
(TE)  for  lipid  suppression  can  cause  artifacts  resulting 
from  subject  motion  and  instrumental  instabilities.  MQ 
filters  alleviate  the  problem  of  motion  artifacts  by  the 
technique  of  a  single-shot  acquisition.  However,  these  se¬ 
quences  generally  exhibit  signal  contamination  from  resid¬ 
ual  lipid  MQ  coherences  (25,26).  In  contrast,  the  selective 
MQ-coherence  (SelMQC)  transfer  sequence  avoids  excita¬ 
tion  of  lipid  MQ  coherences,  and  thus  offers  complete 
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FIG.  1.  The  SS-SelMQC  pulse 
sequence.  The  ZQ->DQ  coher¬ 
ence  transfer  pathway  is  selected 
with  the  Gsei  gradients  in  a  ratio  of 
9i-  02:  g3  =  0:  —  1:  2  (0:  —14 
G/cm:  28  G/cm)  with  duration 
8-1  =  82  =  63  =  2  ms.  All  other 
parameters  are  given  in  Materials 
and  Methods. 


lipid  and  water  suppression  in  a  single  scan  (28).  Like 
most  other  MQ  filters,  the  SelMQC  sequence  has  the  dis¬ 
advantage  of  recovering  only  50%  of  the  lactate  signal,  due 
to  the  selection  of  a  single  coherence-transfer  pathway. 
Although  full  Lac  signal  recovery  is  possible,  a  two-step 
phase  cycling  procedure  is  required  to  refocus  the  lipid 
single  quantum  coherences  and  makes  this  method  sus¬ 
ceptible  to  motion  artifacts.  The  conversion  of  this  edited 
signal  into  accurate  tissue  concentrations  of  Lac  requires 
corrections  for  signal  loss  arising  from  various  additional 
sources.  These  include  molecular  diffusion,  scalar  cou¬ 
pling  (J)  effects,  RF  pulse  imperfections,  spin-lattice  (T}), 
and  spin-spin  (T2)  relaxation  times. 

In  previous  work  (28),  the  SelMQC  technique  has  been 
developed  for  selective  detection  of  Lac  CH3  protons  with¬ 
out  any  interference  from  water  and  co-resonating  lipid 
resonances.  This  technique  has  been  successfully  applied 
for  lactate  detection  in  breast  cancer  (21),  the  detection  of 
an  antineoplastic  agent  Iproplatin  in  murine  RIF-1  tumors 
(27),  and  for  detecting  early  response  to  cyclophospha¬ 
mide  treatment  of  RIF-1  tumors  (17).  Although  this  se¬ 
quence’s  efficiency  for  Lac  editing  with  Lip  suppression 
depends  on  selective  excitation  with  RF  pulses,  it  was 
shown  that  substituting  a  slice  selective  pulse  for  the  first 
frequency  selective  RF  pulse  results  in  tolerable  residual 
lipid  resonances  (28).  To  ensure  complete  lipid  suppres¬ 
sion,  two-step  phase  cycling  of  the  first  CH  pulse  with 
respect  to  the  receiver  is  needed,  which  again  may  have 
motional  instabilities. 

Here  we  report  a  novel  modification  of  the  SelMQC, 
spectral-selective  SelMQC  (SS-SelMQC)  sequence  (Fig.  1) 
using  composite  binomial  pulses  (29)  for  Lac  detection  in 
nonlocalized  and  slice-localization  versions,  with  com¬ 
plete  removal  of  lipid  and  water  resonances.  Lactate  spec¬ 
tra  from  nonlocalized,  one-dimensional  (ID)  localized  sag¬ 
ittal  slice  and  two-dimensional  (2D)  CSI-lactate  images 
were  obtained  from  30  mM  Lac/Lip  phantoms  and  in  vivo 
Dunning  R3327  prostate  animal  tumors  using  our  se¬ 
quence  and  compared  with  the  SelMQC  sequence.  Theo¬ 
retical  simulations  were  generated  to  understand  the  per¬ 


formance  of  the  new  sequence  for  different  parameters 
employed  in  our  experiments. 


In  the  literature  (30),  composite  binomial  pulses  with  in¬ 
creasing  order  were  described  as  {1-1},  (1-2-1),  (1-3-3-1), 
etc.  Initially  binomial  pulses  were  designed  for  spectral- 
selective  excitation  or  inversion  of  spin  resonances  with 
water  suppression  in  NMR  spectroscopy  (30).  In  recent 
years  the  behavior  of  these  pulses  was  explored  for  appli¬ 
cation  to  spatial  selection  in  addition  to  spectral  selection, 
particularly  for  water-fat  MR  imaging  (31-33).  RF  pulses 
and  slice  select  gradients  for  composite  spatial-spectral 
pulses  were  applied  in  different  ways  (31). 

In  our  present  study  we  used  a  (1-1)  composite  pulse 
consisting  of  two  RF  pulses  of  equal  amplitude  with  spac¬ 
ing  (A)  between  them.  Frequency  selective  excitation  was 
achieved  by  using  a  binomial  pulse  block  [(tt/4)+1  —  A 
—  ('rr/4)^>2]  and  adjusting  the  phases  of  RF  pulses  to  selec¬ 
tively  excite  either  off-resonant  Lac  CH3  (spin  I)  reso¬ 
nances  (<j>j,  4 >2  —  x,  — x)  or  on-resonant  CH  (spin  S)  reso¬ 
nances  (<(>!,  <j>2  =  x,  x).  Frequency  selective  inversion  was 
achieved  using  [(tt/  2)x  —  A  —  (ir/  2)  _  J  for  the  Lac  CH3 
resonances  without  perturbing  CH  chemical  shift.  The 
time  difference  between  the  pulses,  “A,”  is  equal  to  the 
reciprocal  of  twice  the  difference  of  the  center  frequencies 
of  maximum  (uj)  and  null  (us)  excitation  bands  of  I  and  S  AQ:  2 
resonances, 

1 

28vJS’ 

where  8u/s  is  the  frequency  difference  I  and  S  resonances. 

When  RF  pulses  are  considered  as  delta  pulses,  A  is  the 
time  delay  between  two  elements  of  the  (1-1)  composite 
pulse.  Under  our  experimental  conditions  where  the  RF 
pulses  are  shaped  pulses  with  finite  pulse  widths,  A  needs 
to  be  corrected  for  finite  pulse  width  effects.  If  we  assume 
Pw(tt/4)  and  pW(„/2j  are  pulse  widths  of  (tt/4)  and  (-tt/2)  RF 
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pulses  within  binomial  90°  and  180°  composite  pulses, 
then  A,  (A  —  pw(lr/4))  and  A2  (A  —  pw(lr/2)l  are  the  effective 
interpulse  delays  (pulse  width  corrected)  for  90°  and  180° 
binomial  blocks,  respectively. 

Numerical  computer  simulations  were  generated  to  cal¬ 
culate  the  selectivity  of  these  pulses  for  the  Lac  “IS”  spin 
system  using  spin  density  matrix  calculations  assuming 
the  experimental  held  strength  of  4.7T.  The  response  of  a 
nuclear  spin  ensemble  to  multiple  RF  irradiation  pulses 
often  deviates  from  the  expected  or  the  desired  effect, 
primarily  due  to  imperfections  such  as  the  spatial  inho¬ 
mogeneity  of  the  RF  field  and  off-resonance  effects.  Due  to 
the  long  durations  of  the  shaped  RF  pulses,  evolution  due 
to  scalar  couplings  and  chemical  shifts  of  both  spins  needs 
to  be  considered.  The  Hamiltonian  of  an  IS  system  during 
the  application  of  a  RF  pulse  and  free  evolution  (no  RF 
effect)  in  the  doubly  rotating  frame  are  given  by 

HRF(t)  =  co1Xt)(4cos4>(t)  +  Lysin4>(f)) 

+  Gols(f)(Sxcos4>(f)  +  Sysin4>(f))  [1] 

and 

H  =  Ao)/z  +  Ad)sSz  +  2rnJISIzSz  [2] 

where  Aooj,  Acos,  JIS  are  chemical  shifts  of  I,  S  spins,  and 
scalar  coupling,  loj/f),  a>is(t ),  and  <\>(t)  are  the  RF  field 
strengths  in  3-lobe  sine  shape  with  step  duration  “t”  on  I 
and  S  spins,  and  their  phase,  respectively.  We  examined 
the  transient  response  of  this  spin  system  following  the 
SS-SelMQC  pulse  sequence  in  the  absence  of  the  slice 
gradient  which  can  be  described  as: 

pitl)  —  1/ total  X  p(0)  X  Utotai  [3] 

where  p(0)  represents  the  equilibrium  density  matrix  and 
the  RF  pulse  effect  was  calculated  from  the  Hamiltonian 
using  Eq.  [1]  by  considering  the  scalar  evolution  term  in 
addition  to  RF  terms  to  account  for  coupling  evolution 
during  finite  pulse  durations.  RF  pulses  within  binomial 
pulse  blocks  individually  are  nonselective  pulses  and 
achieve  spectral  selectivity  as  an  effect  of  composite 
blocks. 

If  we  consider  R(0,4>)  as  a  rotation  propagator  represent¬ 
ing  shaped  RF  pulse  of  flip  angle  0  and  phase  4>  and  U(t)  as 
a  free  evolution  propagator  during  time  period  “t,”  Utotai 
can  be  written  as: 

UtotaI=  U(t2)R[^,x)  UiAJR^x)  t/(t1/2)fl(|,  -  x) 

X  I/(A2)fl(|,xJ  U(tj2)R^  J,x)  E/(A r)J?(  Jx) 

X  U(t)r{^,  -  xj  I7(A,)fl(  J,xj  [4] 

where 

U(t)  =  e~mt  and  R(0,< |>)  =  {hmw  [5] 
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Although  the  propagator  only  deals  with  the  first  two  lines 
of  the  pulse  sequence,  the  effect  of  MQ  gradients  were 
considered  in  the  simulation  by  manually  selecting  the 
MQ  coherence  elements  in  the  density  matrix  to  choose 
the  ZQ^DQ  pathway. 

Hence,  from  Eq.  [3]  the  final  transverse  magnetization  is 
calculated  as: 

<  M^t)  >  =  [  <  Aif(f)  >  +  <  M2( t)  >  ]1/2,  [6] 

where  the  normalized  in-phase  (Mx  (t))  and  out-of  phase 
(My  ( t) )  components  of  the  induced  signal  are  given  as: 

<  Mx(t)  >  =  Tr[p(f2)4]/7rK],  [7] 

<  My(t)  >  =  Tr[p(t2)Iy\/Tr[I^\  [8] 

This  time-domain  complex  signal  is  Fourier  transformed 
to  obtain  the  lactate  spectra  to  test  the  performance  of 
resonance  offsets,  coupling  effects,  and  RF  inhomogeneity. 
We  employed  the  MatLab  software  (MathWorks,  Natick, 
MA)  to  develop  a  home-built  pulse  program  for  simulation 
of  spectral-selective  behavior,  slice  profiles,  to  study  the 
influence  of  scalar  coupling  effects,  RF  finite  pulse  width 
effects,  RF  inhomogeneity  effects,  and  molecular  diffusion 
effects. 

Experimental 

All  MR  imaging  and  spectroscopy  experiments  were  per¬ 
formed  on  a  4.7T  Bruker  Biospin  spectrometer  (40  cm 
horizontal  bore).  Animal  studies  were  conducted  in  com¬ 
pliance  with  protocols  approved  by  the  animal  care  pro¬ 
tocols  in  Memorial  Sloan-Kettering  Cancer  Center. 

Phantom  Preparation 

A  two-compartment  cylindrical  phantom  (25  mm  diame¬ 
ter)  was  prepared  with  30  mM  lactate/H20  doped  with  25 
p,M  Gd-DPTA  and  vegetable  shortening  (Crisco)  side  by 
side.  This  phantom  was  used  to  demonstrate  the  nonlocal- 
ized  and  localized  Lac  editing  with  water  and  lipid  sup¬ 
pression  (Fig.  3).  A  second  phantom  of  10  mM  Lac/D20 
was  prepared  for  generation  of  slice  profile  (Fig.  4)  and 
checking  the  experimental  signal  loss. 

Animal  Preparation  and  Tumor  Volume  Measurements 

Frozen  Dunning  R3327-AT  prostate  cancer  cells  syngeneic 
to  rats  were  thawed  and  seeded  in  75  cm2  flasks.  Copen¬ 
hagen  rats,  supplied  by  Charles  River  Laboratories  (Wil¬ 
mington,  MA),  were  implanted  with  five  million  cells  sub¬ 
cutaneously  in  the  thigh.  Tumor  growth  was  evident  visu¬ 
ally  about  10  days  after  implantation  on  the  thigh  region 
and  volume  increased  to  800  mm3  to  2600  mm3  (not  cor¬ 
rected  for  skin  thickness)  in  an  additional  5-8  days.  Tu¬ 
mor  volume  was  calculated  as  V  =  (tt/6)  *  x  *  y*  z;  where 
x,  y,  and  z  are  the  length,  breadth,  and  depth  of  the  tumor. 

MR  experiments  were  performed  at  a  tumor  volume  of 
857  mm3.  Each  rat  was  anesthetized  using  a  mixture  of 
isoflurane  (1.5-2. 5%)  and  air  and  placed  in  the  animal 
holder.  The  tumor  was  adjusted  to  be  placed  inside  a 
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2-turn  home-built  coil  (25  mm  diameter)  tuned  to  200  MHz 
and  matched  to  50  ohms  prior  to  initiating  MR  experi¬ 
ments.  The  animal  was  placed  inside  a  custom-designed 
MR  platform  with  a  315  mm  long  section  of  120  mm 
diameter  plastic  tube  (the  center  of  this  platform  matches 
with  isocenter  of  the  magnet).  The  animal  was  positioned 
in  the  center  of  this  plastic  cradle  containing  an  RF  coil 
and  was  held  in  place  with  tape.  Temperature  was  main¬ 
tained  at  37°C.  The  magnet  was  shimmed  and  resulted  in 
a  full  width  at  half  maximum  peak  height  (FWHM)  of  less 
than  50  Hz  for  the  water  resonance  in  vivo. 

MR  Imaging 

The  Bruker  ParaVision  Tripilot  pulse  sequence  was  used 
to  create  three  perpendicular  images  (field  of  view  [FOV] 
=  40  mm)  as  scout  images  for  positioning  of  the  animal 
tumor  within  the  coil  to  ensure  that  the  tumor  was  in  the 
center  of  the  magnet.  Using  these  Tripilot  scout  images, 
T2-weighted  MR  sagittal  images  with  a  5  mm  thickness, 
using  a  multislice  multiecho  (MSME)  sequence  (repetition 
time  [TR]  =  3734  ms;  echo  time  [TE]  =  30  ms;  512  X  128 
matrix,  number  of  excitations  [NEX]  =  2;  FOV  =  40  mm) 
were  obtained  from  the  tumor  located  at  isocenter.  Total 
MR  image  acquisition  including  shimming  took  about 
10  min. 

MR  Spectroscopy 

Nonlocalized  Lac  MR  spectra  from  whole  phantom/tu¬ 
mors,  ID  localized  5  mm  sagittal  slice,  as  well  as  2D  CSI 
images  from  a  sagittal  slice  were  collected  using  SS- 
SelMQC  sequence  and  compared  with  the  previously  re¬ 
ported  SelMQC  sequence  (28).  In  vitro  Lac/Lip  phantom 
and  in  vivo  R3327-AT  prostate  tumors  were  studied. 

The  MR  spectroscopy  acquisition  parameters  for 
SelMQC  (28)  included  a  1  ms  three-lobe  sine  shaped  slice- 
selective  pulse  and  a  frequency  selective  15  ms  single-lobe 
Sine  RF  pulses  for  frequency  selective  excitation  and  in¬ 
version  of  CH  and  CH3  resonances.  In  SS-SelMQC  we  used 
a  high-power  three-lobe  sine  shaped  RF  pulses  for  both 
(tt/4)  (200  |jls  duration)  and  (tt/2)  (400  qs  duration)  flip 
angles,  a  pulse  repetition  time  (TR)  of  2  sec,  and  spectral 
width  of  2500  Hz;  t3  =  2  ms;  =  693  qs;  A2  =  493  p,s  (see 
Fig.  1).  Due  to  RF  finite  pulse  widths  within  the  binomial 
pulse  blocks,  A3  and  A2  values  were  adjusted  to  include  the 
chemical  shift  evolution  starting  from  the  center  of  the  first 
pulse  to  the  center  of  the  second  pulse.  In  the  binomial 
spectral-selective  pulses  we  chose  to  use  broadband  sine 
shaped  RF  pulses.  A  slice  gradient  (trapezoidal  shape)  was 
applied  to  choose  a  5  mm  or  10  mm  slice  thickness.  Gra¬ 
dient  rise  time  and  ramp  down  times  were  typically  100- 
150  |jls.  A  phase  cycling  gradient  scheme  with  gp  g2:  g3  = 
0:  —1:  2  (0:  —14  G/cm:  28  G/cm)  with  duration  83  =  S2  = 
S3  =  2  ms  was  employed  in  both  sequences  for  the  selec¬ 
tion  of  ZQ— »DQ  coherence  selection  pathway.  The  trans¬ 
mitter  was  set  at  the  CH  frequency.  The  Lac  CH  resonance 
frequency  is  at  4.1  ppm  compared  to  the  water  signal 
(4.8  ppm).  Whole  tumor  and  ID  slice  localized  Lac  spectra 
were  obtained  with  16  transients,  respectively,  without 
and  with  the  application  of  slice-selection,  and  no  phase 
encoding  gradients  were  applied.  The  2D  3H  CSI  data  were 
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collected  by  applying  16  X  16  phase  encoding  steps  and 
FOV  =  40  mm,  which  results  in  a  voxel  size  of  2.5  X  2.5  X 
5  mm3.  The  total  time  for  the  2D  CSI  experiment  was  =1  hr 
15  min  (NEX  =  8).  The  anatomical  sagittal  T'2-1H  reference 
images  were  obtained  with  the  same  setup  as  the  2D  CSI 
before  starting  the  CSI  scan,  so  the  MR  images  can  be 
coregistered  with  2D  CSI  data.  The  slice  thickness  in  the 
T2-weighted  image  was  matched  with  the  slice  thickness 
in  the  2D-CSI  for  coregistration. 

Data  Processing 

Nonlocalized  and  ID  slice-localized  Lac  spectral  data  were 
processed  using  Bruker  Xwinnmr  software  and  lactate  sig¬ 
nals  were  integrated  using  the  area  under  the  peak.  To 
calculate  the  Lac  signal  enhancement  for  the  SS-SelMQC 
sequence,  we  calculated  the  ratio  of  the  lactate  peak  area  in 
the  SS-SelMQC  to  the  Lac  peak  area  in  the  SelMQC.  Al¬ 
though  quantification  of  lactate  can  be  performed  by  the 
“substitution”  method  (34)  which  was  used  previously 
(35),  since  these  studies  were  done  sequentially  without 
moving  the  sample  (phantom/tumor),  there  was  identical 
coil  loading  and  similar  RF  inhomogeneity  and  therefore 
these  imperfections  will  not  affect  the  Lac  signal  enhance¬ 
ment  ratio.  Signal  enhancement  factors  are  represented  in 
terms  of  mean  ±  standard  deviation.  2D  1H-CSI  images 
were  reconstructed  and  overlaid  with  corresponding  T2- 
weighted  MR  image  using  3DiCSI  processing  software 
(Provided  by  Dr.  Truman  Brown,  Columbia  University).  In 
2D  CSI  data,  voxels  fully  filled  with  Lac/Lip  solution 
(phantom)  were  identified  and  the  lactate  peak  areas 
within  these  voxels  were  calculated  using  iterative  non¬ 
linear  least-squares  method  using  XSOS  package  (pro¬ 
vided  by  Dr.  Dikoma  Shungu,  Weill  Medical  College).  The 
exact  same  voxels  were  analyzed  from  the  data  obtained 
using  both  sequences  without  considering  inhomoge¬ 
neity  factors  and  coil  loading  factor. 

RESULTS 

Excitation  and  inversion  profiles  following  binomial 
pulses  were  generated  using  an  in-house  computer  pro¬ 
gram  representing  the  experimental  conditions.  The  selec¬ 
tive  excitation  profile  for  the  binomial  spectral-selective 
pulse  is  shown  for  excitation  of  —2.8  ppm  resonance  with¬ 
out  disturbing  0  ppm  resonance  (Fig.  2a)  and  vice  versa  F2 
(Fig.  2b).  Similarly,  the  selective  inversion  profile  of  the 
binomial  spectral-selective  pulse  excites  the  “  —  2.8”  ppm 
resonance  without  disturbing  “0”  ppm  resonance  (Fig.  2c). 

In  generating  this  figure,  the  Lac  CH  peak  was  referenced 
to  “0”  ppm  and  hence  Lac  CH3  will  appear  at  “2.8”  ppm  on 
the  chemical  shift.  Figure  3  a  shows  the  MR  image  of  a  F3 
5-mm  thick  coronal  slice  across  the  phantom  showing  the 
two  compartments  with  Lac  and  Lip,  respectively.  Figure 
3b  shows  the  nonlocalized  Lac  spectra  obtained  using 
SS-SelMQC  and  compared  with  SelMQC  (Fig.  3c)  from  the 
whole  phantom  with  excellent  lipid  suppression.  Lipid 
suppression  was  verified  by  choosing  a  slice,  using  SS- 
SelMQC  and  SelMQC  from  the  lipid  portion  of  the  cylin¬ 
drical  phantom,  and  no  signal  was  detected  at  the  Lac/Lip 
frequency  (data  not  shown  as  voxels  spectra  contain  no 
signal).  Lac  signal  enhancement  was  observed  in  SS- 
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FIG.  2.  Theoretical  excitation  following  |‘(n/4)+1  —  \  —  (ir/4)^2) 
pulse  sequence  for  spectral  selection  with  (4>j,  ij>2  =  x,  -x)  (a)  and 
(4>i,  4>2  =  x,  x)  (b)  as  a  function  of  resonance  offset.  Inversion  profile 
of  for  inversion  of  -2.8  ppm  resonances  shown  in  (c).  In  these 
simulations  the  Lac  CH  peak  was  referenced  to  “0”  ppm  and  hence 
Lac  CH3  will  appear  at  “2.8”  ppm  on  the  chemical  shift.  RF  pulse 
parameters  and  interpulse  durations  were  chosen  to  match  the 
experimental  conditions  given  in  the  text. 


FIG.  4.  Experimental  spatial-profile  of  lactate  following  the  slice 
selective  RF  pulse  block.  TR  =  3  sec,  number  of  scans  =  1 ,  number 
of  data  points  =  512,  total  acquisition  time  less  than  5  min.  We  used 
10  mM  Lac/D20  phantom  for  this  experiment. 


SelMQC  over  SelMQC  with  no  significant  difference  in 
performance  with  respect  to  water  and  Lip  suppression 
levels.  We  repeated  this  experiment  at  three  different 
times  (each  timepoint  corresponds  to  a  new  experiment 
starting  with  positioning,  tuning,  and  matching)  to  test  the 
consistency  of  signal  enhancement  in  lactate  signal  from 
the  whole  Lac/Lip  phantom.  The  integrated  lactate  signal 
area  in  SS-SelMQC  is  2  to  3  times  the  lactate  signal  ob¬ 
served  in  SelMQC  with  similar  suppression  of  lipid  and 
water  signals.  Using  a  paired  f-test  from  EXCEL  software, 
the  differences  in  the  Lac  signal  were  statistically  signifi¬ 
cant  ( P  <  0.05)  between  different  experiments  with  the 
signal  enhancement  factor  being  2.4  ±  0.40  in  nonlocal- 
ized  SS-SelMQC,  compared  to  SelMQC.  Similarly  we 
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tested  localized  in  vitro  ID  Lac  signal,  which  yielded 
similar  signal  enhancement  factors  (data  not  shown).  Lo¬ 
calization  was  achieved  by  applying  slice-select  gradient 
on  the  first  spectral-selective  CH3  binomial  block  and  its 
spatial  profile  is  shown  in  Fig.  4. 

The  efficacy  of  this  pulse  sequence  was  also  analyzed  by 
generating  lactate  signals  derived  from  localized  slices.  2D 
chemical  shift  imaging  data  were  obtained  from  a  10-mm 
sagittal  slice  across  the  phantom  located  at  the  center 
(covering  both  Lip  and  Lac)  using  both  methods  (data  not 
shown).  With  similar  water  and  lipid  suppression,  lactate 
signal  from  different  voxels  shows  a  statistically  signifi¬ 
cant  (P  <  0.05)  signal  enhancement  factor  of  2.4  ±  0.29 
times  in  SS-SelMQC,  compared  with  the  SelMQC  se¬ 
quence.  The  result  is  similar  to  the  enhancement  observed 
in  the  nonlocalized  spectra  obtained  using  both  sequences. 

The  same  effect  was  demonstrated  using  R3337  tumors. 
Figure  5a  shows  the  stacked  spectrum  of  lactate  signal 
obtained  from  a  5-mm  sagittal  slice  from  a  tumor  (vol¬ 
ume  =  857  mm3)  using  both  sequences.  From  peak  area 
integrals,  the  lactate  signal  to  noise  with  SS-SelMQC  is 
more  than  200%  that  obtained  with  the  SelMQC.  We  also 
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FIG.  3.  a:  Experimental  coronal  MR  image  obtained  from  lactate/ 
lipid  two-compartment  phantom.  Unlocalized  proton  spectra  of  lac¬ 
tate  obtained  from  this  phantom  were  compared  using  SS-SelMQC 
(b)  and  SelMQC  (c).  Enhancement  of  2.8  times  was  observed  in 
SS-SelMQC  compared  with  SelMQC.  Water  suppression  with  a  new 
sequence  is  similar  to  that  obtained  with  SelMQC.  Data  were  ac¬ 
quired  with  1 6  scans  with  TR  =  2  sec.  spectral  width  =  2500  Hz  and 
total  acquisition  time  was  less  than  2  min.  Measured  linewidths  of 
Lac  signal  were  34  Hz,  29  Hz  in  (b)  and  (c),  respectively. 


FIG.  5.  a:  Stacked  plot  of  experimental  ID  localized  %  lactate 
spectrum  obtained  from  in  vivo  R3327-AT  tumor  (volume  = 
857  mm3)  using  SS-SelMQC  and  SelMQC.  Experimental  parame¬ 
ters  are  described  in  the  text,  b:  2D  Lactate  CSI  from  sagittal  5-mm 
thick  slice  using  SS-SelMQC.  CSI  data  is  overlaid  on  T2-contrast 
MR  image.  Chemical  shift  scale  displayed  over  the  tumor  is  between 
2. 5-0. 5  ppm.  Measured  linewidths  of  Lac  signal  were  27  Hz,  25  Hz 
using  SelMQC  and  SS-SelMQC  sequences.  Number  of  scans  =  8, 
TR  =  2  sec,  and  total  acquisition  time  1  hr  15  min. 
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obtained  16  X  16  2D-CSI  lactate  maps  (Fig.  5b)  coregis¬ 
tered  with  a  T2-weighted  image  corresponding  to  the  slice 
shown  in  Fig.  5a.  For  comparison  of  both  sequences,  we 
did  not  move  the  sample  between  the  two  experiments. 

DISCUSSION 

Enhanced  signal-to-noise  ratios  of  metabolites  are  impor¬ 
tant  since  it  will  lead  to  decreased  scan  time,  which  is 
important  both  in  clinical  studies  to  minimize  motion  and 
patient  discomfort,  and  in  preclinical  studies  to  minimize 
anesthesia,  or  higher  spatial  resolution. 

Comparison  of  SS-SelMQC  and  SelMQC  Pulse 
Sequences 

Following  a  series  of  preparation  RF  pulses  and  J-evolu- 
tion  delays  in  both  sequences,  lipid  protons  in  addition  to 
lactate  protons  can  generate  MQ  coherences.  By  applying 
Lac  CH-  and  CH3-  frequency  selective  excitation  pulses 
and  CF13-  inversion  pulse,  the  fraction  of  lipid  MQ  coher¬ 
ences  that  can  be  converted  to  observable  magnetization 
was  minimized.  In  SelMQC,  frequency  selective  sine 
pulses  were  applied,  whereas  in  SS-SelMQC  these  pulses 
were  replaced  by  binomial  {1-1}  pulses.  We  chose  the 
ZQ^DQ  pathway  in  both  sequences  resulting  in  observa¬ 
tion  of  only  50%  Lac  signal.  Although  full  Lac  signal 
recovery  is  possible  by  adding  ZQ— »DQ  and  DQ^ZQ  co¬ 
herence  transfer  pathways  using  a  two-step  phase  cycling 
of  Lac  CH  pulse  (28),  this  modification  would  make  the 
method  susceptible  to  subject  motion  and  instrumental 
instabilities. 

Single-slice  localization  in  the  SelMQC  sequence  by  re¬ 
placing  the  first  frequency  selective  CH3-  excitation  pulse 
with  a  broadband  excitation  3-lobe  sine  pulse  may  intro¬ 
duce  residual  MQ  coherences,  although  it  was  noted  by  He 
et  al.  (28)  that  this  degradation  in  lactate  editing  efficiency 
is  not  noticeable.  This  effect  was  avoided  in  SS-SelMQC 
by  replacement  of  the  slice  selective  RF  pulse  of  SelMQC 
with  a  binomial  pulse  that  removes  the  generation  of  re¬ 
sidual  lipid  single  quantum  coherences  that  interferes 
with  lactate  signal  and  therefore  reduces  the  lactate  signal 
contamination.  The  theoretical  profile  for  spectral  selec¬ 
tion  (based  on  Fig.  2)  has  sharp  transitions  with  «=50  Hz 
range  and  is  within  the  nonexcitation  band.  For  our  study, 
prostate  tumors  were  implanted  on  the  thigh  region  and 
respiration  implied  frequency  shift  is  likely  to  be  negligi¬ 
ble.  But  we  propose  to  use  higher-order  binomial  pulses  to 
have  a  broader  nonexcitation  region  to  alleviate  the  spec¬ 
tral  selectivity  issues  due  to  the  frequency  shifts  that  occur 
due  to  respiration.  In  SS-SelMQC,  by  manipulating  the 
slice  gradient  profile,  this  binomial  pulse  acts  as  a  spectral- 
spatial  pulse,  although  higher-order  binomial  pulses  are 
desirable  for  improved  spectral-spatial  behavior. 

Reasons  for  Signal  Enhancement  in  SS-SelMQC 

Reduced  effects  of  J-scalar  coupling  evolution,  molecular 
diffusion,  as  well  as  T2  relaxation  loss  in  SS-SelMQC  are 
the  likely  causes  of  the  signal  enhancement  observed  ex¬ 
perimentally  compared  to  SelMQC.  These  effects  were 
simulated  using  computer  simulations  of  unlocalized  lac- 
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tate  signals  assuming  an  IS  spin  system.  Lac  editing  effi¬ 
ciency  in  SelMQC  sequence  depends  on  frequency  selec¬ 
tive  excitation  or  inversion  profiles,  which  necessitates  the 
use  of  long  pulses  which  increases  the  MQ  evolution  pe¬ 
riod  and  affects  the  lactate  editing  efficiency  (28).  In  SS- 
SelMQC  we  replaced  all  the  long  frequency-selective 
pulses  with  short  binomial  composite  RF  pulses  for  fre¬ 
quency  selection,  which  resulted  in  shorter  MQ  evolution 
delays.  This  effect  increased  the  Lac  editing  efficiency  by 
reducing  J-coupling  evolution  as  well  as  molecular  diffu¬ 
sion  effects  as  described  below. 

Scalar  Coupling  Effects 

Approximately  24%  of  the  increase  in  the  lactate  signal 
observed  in  SS-SelMQC  (Fig.  6a)  in  comparison  with  F6 
SelMQC  (Fig.  6b)  (assuming  no  RF  pulse  imperfections) 
can  be  explained  from  increased  scalar  coupling  evolution 
during  RF  pulse  finite  widths.  We  denoted  RF  pulse  im¬ 
perfection  (Wj/W?),  where  and  Wj  are  actual  and  ideal 
RF  field  strength  experienced  by  samples  within  RF  coils. 
Hence  (Wj/M/?)  =  1  corresponds  to  no  RF  pulse  imperfec¬ 
tions  and  (W]/W?)  =  0.9,  0.8  corresponds  to  the  presence  of 
10%  and  20%  RF  pulse  imperfections,  respectively.  In  Fig. 

6a, b  the  X-axis  is  the  chemical  shift  scale  and  the  lactate 
doublet  is  centered  on  the  reference  chemical  shift  “1.32” 
ppm,  which  is  the  experimental  frequency  of  CH3.  The 
peak  intensity  is  plotted  with  (J  0)  and  without  (J  =  0) 
scalar  coupling  evolution  considered  during  RF  finite 
pulse  widths  in  the  presence  of  RF  pulse  imperfections. 


FIG.  6.  Computer  simulated  lactate  spectra  generated  as  a  function 
of  RF  pulse  imperfections  and  scalar  coupling  evolutions  (J  Y  0) 
during  RF  pulses  using  (a)  SS-SelMQC  and  (b)  SelMQC  sequence. 
Parameters  used  for  generating  these  line  plots  are  marked  on  the 
figures.  Y-axis  represents  percentage  lactate  signal  normalized  to 
100%  when  there  were  no  scalar  coupling  effects  (J  =  0)  or  RF 
imperfections.  These  graphs  show  simulated  spectra  of  lactate  CH3 
resonance  corresponding  to  the  4.7T  experimental  conditions  and 
assuming  IS  spin  system  for  lactate  (CH  =  I  and  CH3  =  S)  with 
scalar  coupling  J,s  =  7  Hz. 
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The  Y-axis  denotes  the  Fourier  transform  of  total  trans¬ 
verse  magnetization  “Mxy”  following  the  SS-SelMQC  and 
SelMQC  sequences.  Lac  signal  intensity  was  normalized  as 
100%  when  there  were  no  scalar  coupling  effects  and  RF 
imperfections.  Introducing  scalar  evolution  during  finite 
pulse  widths  causes  a  decrease  in  Lac  signal  intensity  to 
99%  in  SS-SelMQC,  with  further  reductions  to  87%  with 
10%  and  59%  with  20%  RF  pulse  imperfections.  Simi¬ 
larly,  in  SelMQC,  Lac  signal  intensity  decreases  to  75%, 
further  reduces  to  62%  with  10%  and  40%  with  20%  RF 
pulse  imperfections.  These  simulations  yielded  a  24%  sig¬ 
nal  gain  for  SS-SelMQC  compared  with  SS-SelMQC.  It 
should  be  noted  that  for  patient  studies  at  the  clinical 
magnetic  field  strengths  of  1.5T  or  3.0T  the  signal  loss  due 
to  scalar  coupling  effects  will  be  more  severe  using 
SelMQC  compared  to  SS-SelMQC.  The  reason  for  this  is 
that  frequency  selective  pulses  used  in  SelMQC  need  to  be 
very  long,  as  chemical  shift  dispersion  is  very  low  and 
results  in  longer  ti  MQ  evolution  delays  leading  to  higher 
signal  loss.  Relatively  much  shorter  L  values  can  be  main¬ 
tained  by  binomial  selective  pulses  in  the  SS-SelMQC 
sequence  resulting  in  minimum  scalar  coupling  losses  in¬ 
creases  Lac  signal  sensitivity.  We  can  further  minimize 
signal  loss  observed  in  simulations  due  to  RF  imperfec¬ 
tions  by  customized  RF  coils  generating  homogeneous  RF 
fields. 

T2  Relaxation  Losses 

Relatively  short  T2  values  of  lactate  have  been  reported  in 
in  vivo,  for  e.g.,  rat  C6-glioma  with  200  ms  (36)  and  mice 
tumors  (37)  with  MCa  (68  ms)  and  Colon-38  tumors 
(117  ms),  which  would  lead  to  lower  Lac  detection  effi¬ 
ciency.  Since  the  in  vivo  T2  relaxation  times  are  short,  an 
optimal  sequence  minimizing  T2  losses  is  critical  for  ac¬ 
curate  measurements  and  for  absolute  quantification  of 
tissue  lactate  concentrations.  This  effect  may  be  particu¬ 
larly  important  in  studies  with  small  tumor  volumes, 
where  lactate  concentrations  maybe  low  (35),  or  if  used  for 
early  cancer  diagnosis.  Since  the  total  duration  of  RF 
pulses  in  the  modified  sequence  are  shorter  and  ==45  ms 
less  than  SelMQC  sequence,  this  enhances  the  lactate  sig¬ 
nal  by  20%  (assuming  T2  =  200  ms)  in  the  SS-SelMQC 
method  by  further  reducing  T2  relaxation  losses. 

Molecular  Diffusion  Effects 

When  two  gradients  are  applied  during  the  MQ  evolution 
period  to  choose  ZQ-^DQ  pathways,  Lac  signal  attenua¬ 
tion  occurs  due  to  molecular  diffusion  of  molecules  during 
the  period  (28)  between  these  gradients  of  duration  8  ms, 
following  the  Stejskal-Tanner  equation  (38).  We  assumed 
phantom  conditions  in  our  calculations  using  the  diffu¬ 
sion  coefficient  D  =  1  X  ICC5  cm2  s-1  at  room  temperature. 
Using  the  equation  reported  earlier  (28),  the  signal  loss  due 
to  molecular  diffusion  is  as  follows: 


S/S(  0)  =  exp 


D(p(  2/TTySg)}2 


[9] 


Where  A  =  L  +  8  +de,  p  is  the  coherence  order,  de  is 
gradient  dead  time,  y  is  gyromagnetic  ratio  =  2.675e8 


s- 1T_1;  for  our  present  SelMQC  parameters,  the  lactate 
signal  loss  will  be  30%  and  signal  loss  in  SS-SelMQC  will 
be  only  6%.  Hence  a  signal  gain  of  24%  is  obtained  in 
SS-SelMQC.  The  molecular  diffusion  losses  will  be  severe 
in  choosing  the  DQ^-ZQ  pathway  to  obtain  the  full  signal 
intensity  of  lactate,  or  at  low  magnetic  field  strengths 
where  long  pulses  result  in  higher  signal  loss. 

Although  the  signal  loss  was  calculated  theoretically  for 
individual  imperfections  by  assuming  other  imperfections 
as  zero,  much  higher  signal  loss  obtained  in  our  experi¬ 
ments  may  be  due  to  these  combined  imperfection  effects 
simultaneously  in  addition  to  magnetic  Held  inhomogene¬ 
ity  effects.  Signal  enhancement  was  also  confirmed  from 
the  Lac  signal  areas  from  SS-SelMQC,  a  single  hard  pulse, 
and  SelMQC  scans.  We  prepared  10  mM  Lac/D20  solution 
as  a  phantom  to  test  the  nonlocalized  lactate  signal  ob¬ 
tained  by  one  pulse  (rectangular  shape),  SelMQC,  and 
SS-SelMQC  sequences.  If  we  normalize  lactate  signal  from 
single  pulse  to  100%,  «=50%  lactate  signal  is  obtained  from 
SS-SelMQC  compared  to  24%  signal  obtained  from 
SelMQC  using  identical  experimental  parameters  for  the 
latter  sequences  (unpubl.  data).  The  possible  in  vivo  lac¬ 
tate  detection  limit  of  SS-SelMQC  sequence  is  in  the  range 
of  2-4  mM. 


CONCLUSIONS 

In  summary,  the  demonstrated  ability  of  the  SS-SelMQC 
sequence  to  measure  tumor  lactate  in  the  presence  of  sub¬ 
cutaneous  lipids  facilitate  quantitative  mapping  of  tumor 
Lac  distributions  in  in  vivo  tumors.  With  lipid  and  water 
completely  suppressed  in  a  single  scan,  the  lactate  signal 
can  be  detected  by  SS-SelMQC  from  a  desired  sample  or 
slice  of  interest  in  tumor  tissues  that  contain  lipid.  These 
results  demonstrate  that  the  SS-SelMQC  sequence  can  de¬ 
tect  lactate  with  higher  sensitivity  without  compromising 
the  editing  efficiency  of  SelMQC  with  similar  lipid  and 
water  suppression  in  a  single  scan.  Lactate  signal  obtained 
from  whole  tumors  detected  by  SS-SelMQC  has  a  statisti¬ 
cally  significant  (P  <  0.05)  signal  enhancement  of  2.4  ± 
0.40  times  compared  to  the  SelMQC  sequence.  ID  local¬ 
ized  and  2D  CSI  maps  were  obtained  from  5-mm  thick 
sagittal  slice  using  SS-SelMQC  and  showed  similar  sensi¬ 
tivity  enhancement.  These  preliminary  studies  indicate 
efficient  selection  of  a  slice  element  and  clearly  point  out 
the  possibilities  for  obtaining,  via  the  methodology  out¬ 
lined  here,  volume  selection  such  that  the  desired  metab¬ 
olite  region  of  interest  is  excited  without  exciting  the  un¬ 
wanted  lipid  and  water  resonances.  We  are  in  the  process 
of  applying  higher-order  binomial  pulse  such  as  (1-3-3-1) 
for  uniform  spectral-spatial  behavior  for  slice  and  volume 
excitations.  This  method  can  be  useful  for  carrying  out 
studies  for  detecting  increased  glycolysis  in  exercised 
muscles  and  human  breast  cancer,  as  well  as  for  lactate 
detection  in  the  internal  organs  of  transgenic  animal  tumor 
models. 
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29  Abstract 

30  We  performed  longitudinal  lactate  spectroscopic  imaging  and  dynamic  contrast  enhanced 

31  magnetic  resonance  imaging  (DCE-MRI)  studies  in  two  prostate  tumor  models,  the  Dunning 

32  R3327  AT  and  Dunning  H,  to  determine  the  potential  of  lactate  and  DCE-MRI  as  markers  of 

33  tumor  aggressiveness.  Dunning  R3327-AT  (16  rats)  and  Dunning  H  (6  rats)  were  studied  at  five 

34  different  tumor  volumes  (groups  I  to  V),  ranging  in  size  from  110-2453  mm3).  Lactate 

35  determinations  were  perfonned  with  the  Selective  Multiple  Quantum  Coherence  (Sel-MQC) 

36  sequence.  Tumor  perfusion  was  evaluated  using  DCE-MRI.  Lactate  was  not  detected  in  small 

37  (110-265  mm3)  R3327-AT  tumors.  Tumor  lactate  increased  from  1.53  ±  0.54  mM  (Group 

38  11=478-649  mm3)  to  4.78  ±  1.66  mM  (Group  III;  1035-1207  mm3),  4.64  ±  1.26  mM  (Group  IV; 

39  1409-1708  mm3),  and  decreased  to  2.71  ±  0.93  mM  (Group  V;  1941-2453  mm3)  and  was  greater 

40  in  the  core  than  the  rim.  Lactate  was  not  detected  in  Dunning  H  tumors.  Akep  decreased  in  both 

41  tumors  with  volume  (p<0.001),  and  was  greater  in  the  rim  than  the  core  (R3327-AT  (p<0.001); 

42  Dunning  H  (p<0.023)).  Dunning  H  tumors  had  higher  Akep  values  than  R3327-AT  tumors 

43  (p<0.001)  demonstrating  better  perfusion.  Lactate  is  not  detected  in  small  Dunning  R3327  AT 

44  tumors  but  increases  with  growth  and  subsequently  decreases.  This  increase  in  lactate  in  the 

45  Dunning  R3327-AT  occurs  with  concomitant  decreases  in  Akep  and  evolving  differences  in 

46  perfusion  between  the  rim  and  core,  indicating  poorer  perfusion.  Lactate  is  not  detected  in  slow 

47  growing  Dunning  H  tumors  suggesting  it  as  a  potential  biomarker  for  aggressive  tumors. 

48 

49  Abbreviations:  2D-CSI,  2  dimensional-chemical  shift  imaging;  DCE-MRI,  dynamic  contrast- 

50  enhanced  magnetic  resonance  imaging;  HIF,  hypoxia  inducible  factor;  MRSI,  magnetic 

5 1  resonance  spectroscopic  imaging;  Sel-MQC,  selective  multiple  quantum  coherence. 
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52  Introduction 

53  Patients  with  the  identical  pathological  cancer  diagnosis  often  have  different  clinical  courses. 

54  This  unexplained  heterogeneity  is  particularly  true  of  prostate  cancer,  a  disease  which  is  usually 

55  indolent,  yet  sometimes  follows  an  aggressive  course.  Discriminating  between  patients  with 

56  aggressive  vs.  indolent  prostate  cancer  is  a  major  clinical  problem  and  often  leads  to  patients 

57  being  overtreated  [1-2].  Nomograms  are  useful  to  predict  outcomes  and  for  selecting  patients  for 

58  watchful  waiting  [3-4]  but  are  imperfect,  and  newer  biomarkers  with  higher  predictive  value  are 

59  necessary  to  add  to,  and  improve  upon  clinical  parameters  (PSA,  Gleason  grade,  stage  etc). 

60 

61  In  addition  to  tumor  heterogeneity  between  patients,  within  a  specific  tumor  there  are  metabolic, 

62  physiologic  and  genetic  differences  between  cells,  which  can  cause  treatment  failure  and 

63  metastases.  Clonal  evolution  and  changes  in  the  tumor  microenvironment  during  tumor  growth 

64  can  lead  to  the  development  and  selection  of  cells  with  more  malignant  phenotypes  [5-6]. 

65  Changes  in  tumor  physiology  can  induce  hypoxia  which  can  cause  preferential  selection  of  cells 

66  with  p53  mutations  and  greater  metastatic  potential  [7-8],  tumor  progression  [9],  and  resistance 

67  to  anti-neoplastic  treatments  [10-11],  Detecting  local  areas  of  resistance  to  treatment  with  new 

68  more  accurate  biomarkers  would  allow  targeting  of  these  areas  with  local  therapies  or  higher 

69  radiation  doses  [12],  to  eradicate  resistant  cells.  Lactate  has  been  recognized  in  cervical  cancer 

70  and  head  and  neck  tumors  as  a  marker  of  a  higher  risk  of  metastases  and  poor  prognosis  [13-14]. 

71 

72  Changes  in  lactate  have  also  been  suggested  to  be  an  early  marker  of  response  [15],  [16]. 

73  Tumors  evolve  with  growth  and  therefore  a  non-invasive  biomarker  might  be  important  for 

74  repetitive  non-invasive  assays  of  the  tumor.  With  the  widespread  use  of  MRI  and  emerging  use 
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75  of  MRSI  in  prostate  cancer  staging  [17-18],  a  modification  to  existing  pulse  sequences  could 

76  readily  non-invasively  detect  lactate. 

77  The  Dunning  R3327-AT  tumor  is  an  aggressive,  highly  anaplastic  prostate  cancer  with  a 

78  doubling  time  of  3.3  days  which  was  derived  from  a  Dunning  R3327-H  tumor,  a  well- 

79  differentiated,  hormone  responsive  tumor  with  a  doubling  time  of  12-20  days.  The  well- 

80  differentiated,  slower  growing  Dunning  H  has  higher  oxygen  levels  [19]  compared  to  anaplastic 

81  R3327  AT  tumors,  longer  radiation  induced  growth  delays  [20],  and  lower  pimonidazole  binding 

82  (5%)  compared  to  R3327-AT  tumor  (18%;  p<0.05%)  [21].  Therapeutic  and  histological 

83  characteristics  of  these  tumors  are  similar  to  human  prostate  adenocarcinomas  [20]. 

84 

85  Imaging  lactate  in  vivo  by  magnetic  resonance  is  feasible  but  challenging  since  the  methyl 

86  resonance  of  lactate  is  co-resonant  with  the  methyl  resonance  of  lipids.  The  scalar  (J)  coupling 

87  between  the  doublet  at  1.3  ppm  arising  from  the  methyl  (CH3)  protons  and  the  quartet  at  4.1  ppm 

88  arising  from  the  methine  proton  (CH)  is  used  to  selectively  detect  the  lactate  signal  [22-24].  The 

89  selective  multiple-quantum  coherence  transfer  (Sel-MQC)  [25]  sequence  with  phase  encoding 

90  was  used  to  generate  lactate  maps. 

91 

92  The  aim  of  this  study  was  to  compare  the  concentration  of  lactate  in  aggressive  and  indolent 

93  prostate  tumor  model  to  determine  its  potential  as  a  marker  of  tumor  aggressiveness. 

94  Specifically,  the  goals  were  to  (a)  quantitatively  compare  the  whole  tumor  lactate  concentration 

95  in  Dunning  R3327-AT  and  Dunning  H  prostate  tumor  models  to  determine  the  potential  of 

96  lactate  as  a  marker  of  tumor  aggressiveness,  (b)  evaluate  spatial  heterogeneity  of  lactate  and  c) 

97  evaluate  tumor  perfusion  and  penneability  in  both  tumor  models  using  dynamic  contrast 


98 


enhanced  MRI  (DCE-MRI). 
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99  Experimental  Methods 

100  Animal  model  and  cell  culture 

101  Animal  studies  were  conducted  in  compliance  with  protocols  approved  by  the  Institutional 

102  Animal  Care  and  Use  Committee  (IACUC)  of  Memorial  Sloan-Kettering  Cancer  Center 

103  (MSKCC).  Dunning  R3327-AT  cells  were  seeded  in  75  cm2  flask  and  grown  in  20  ml 

104  Dulbecco’s  Modified  Essential  Medium  (DMEM)  containing  4.5  g/L  of  Glucose  and  L- 

105  glutamine  without  sodium  pyruvate  and  supplemented  with  FCS  (10%),  Penicillin  and 

106  Streptomycin  (1  %)  were  added  under  sterile  conditions  at  37  °C  in  a  95  %  humidified  air  and  5 

107  %  carbon  dioxide  atmosphere.  Trypsin  (0.05%)  with  0.53  inM  EDTA,  HBSS  without  calcium 

108  and  magnesium  salts  was  used  to  separate  the  cells.  Copenhagen  rats  (Charles  River 

109  Laboratories)  were  implanted  with  5xl06  cells  subcutaneously  in  the  thigh  region.  Dunning  H 

110  tumor  was  generously  provided  by  Drs.  Susan  Dalrymple  and  John  Isaacs  from  Johns  Hopkins 

111  School  of  Medicine  and  tumors  were  grown  by  implanting  2  mm  sections  of  tissue 

112  subcutaneously  in  the  thigh  of  Fischer-Copenhagen  rats  (Harlan  Sprague  Dawley). 

113 

114  Tumor-bearing  rats  were  studied  longitudinally  and  the  data  grouped  based  on  measurements 

115  performed  at  similar  tumor  volumes.  For  the  Dunning  R3327-AT  tumor,  the  groups  were  1 10- 

1 16  265  mm3  (I),  478-649  mm3  (II),  1035-1207  mm3  (III),  1409-1708  mm3  (IV)  and  1941-2453  mm3 

•5 

1 17  (V).  For  the  Dunning  H  models,  the  studies  were  similarly  grouped  as  227-259  mm  (I),  430-546 

118  mm3  (II),  1000-1236  mm3  (III),  1382-1715  mm3  (IV)  and  1886-2255  mm3  (V).  Tumor  volume 

119  (V)  was  calculated  as  V=  (k/6)*  x  *  y  *  z;  where  x,  y  and  z  are  the  length,  breadth  and  vertical 

120  depth  of  the  tumor,  correcting  for  the  skin  thickness  of  1.3  mm  [20] 

121 
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122  MR  Experiments 

123  Tumor  bearing  rats  were  anesthetized  with  Isoflurane  (1.0  -  2.5  %)  combined  with  air  during 

124  catheterization  and  MR  experiments.  A  tail  vein  catheter  (24G  X  %”  (Terumo  Medical 

125  Corporation,  Elkton,  M.D)  was  inserted  prior  to  placing  the  rat  in  the  magnet.  MR  experiments 

126  were  performed  on  a  Bruker  4.7  T,  40  cm  bore  animal  scanner.  A  home-made  2  turn  solenoid  (25 

127  mm  diameter)  was  used  as  transmitter/receiver  for  all  studies.  Three  perpendicular  scout  images 

128  were  obtained  (Field  of  view  (FOV)=15  cm,  slice  thickness  5  mm,  repeat  interval  (TR)  =  100 

129  ms,  echo  time  (TE)  =  5  ms)  for  positioning  of  the  animal  to  ensure  that  the  tumor  was  in  the 

130  center  of  the  magnet.  Using  the  tri-pilot  scans  as  a  spatial  reference,  a  spin-spin  (T2)  weighted 

131  MR  image  (sagittal  orientation)  of  5  mm  slice  thickness  was  collected  using  the  Multi  Slice 

132  Multi  Echo  (MSME)  sequence  (TR-3734  ms,  TE-  30ms,  a  512  X  128  matrix  size,  1  slice  and  2 

133  excitations). 

134 

135  Lactate  detection  was  perfonned  using  the  Selective  Multiple  Quantum  Coherence  (Sel-MQC) 

136  editing  sequence  [25],  using  frequency  selective  15  ms  single-lobe  Sine  pulses.  The  zero 

137  quantum  (ZQ)— ^double  quantum  (DQ)  coherence  transfer  pathway  was  selected  using  a  phase 

138  cycling  gradient  scheme  with  gl:g2:g3  =  0:-l:2  where  gl,  g2,  and  g3  are  the  relative  gradient 

139  strengths  with  respective  durations  8i  =  82  =  2  ms,  83  =  4  ms,  and  an  amplitude  of  24  G/cm.  Five 

140  hundred  twelve  data  points  were  collected  with  8  averages,  TR=2s,  and  spectral  width  of 

141  2500Hz.  Localization  was  achieved  using  a  16x16  phase  encoding  matrix,  and  a  FOV=40  mm 

142  (2. 5x2. 5  mm  in  plane  resolution).  The  voxel  volume  was  31.2  mm  and  total  2D-CSI  (two 

143  dimensional  chemical  shift  image)  acquisition  time  was  75  minutes.  Two  dimensional  CSI  maps 

144  of  lactate  were  obtained  from  a  5  mm  slice  in  the  tumor  which  was  coregistered  with  the  5  mm 
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3 

145  thick  T2-weighted  image.  For  quantitation  using  the  substitution  method  [26],  a  9.12  cm 

146  cylindrical  reference  phantom  containing  30  mM  lactate  and  25uM  Gd-DTPA  was  scanned  with 

147  identical  parameters. 

148 

149  Upon  completion  of  the  MRSI  lactate  study,  DCE-MRI  imaging  was  done.  A  pre-contrast  proton 

150  density  sagittal  MR  image  (3  slices,  slice  gap=0.20  mm)  was  obtained  by  using  a  Gradient  Echo 

151  Fast  Imaging  (GEFE)  sequence  with  TR/TE=500ms/3ms,  30°  flip  angle,  2  mm  slice  thickness, 

152  4cm  FOV,  128  X  128  matrix  size,  and  2  excitations.  Subsequently  the  dynamic  spin  lattice  (Tl) 

153  weighted  MR  images  were  obtained  using  the  same  sequence/parameters  as  in  the  pre-contrast 

154  proton  density  images  except  that  the  TR  was  reduced  to  50ms,  12  sec  temporal  resolution,  96 

155  time  points,  and  30°  flip  angle.  Gd-DTPA  (0.2  mM/kg;  Magnevist,  Berlex  Laboratories)  was 

156  injected  intravenously  after  the  first  10  images  of  the  dynamic  sequence  were  acquired  and  these 

157  images  were  used  to  measure  the  pre-contrast  signal  intensity  under  saturation  conditions. 

158 

159 

160 

1 6 1  MR  Data  Processing 

162  DCE-MRI  data  were  analyzed  using  the  two  compartment  model  of  Hoffman  et  al  [27]  to 

163  evaluate  the  rate  constant  kep  (rate  constant  for  Gd  moving  from  the  interstitial  space  to  plasma) , 

164  Akep  (A  -  Amplitude  in  the  model)  on  a  voxel-by-voxel  basis  using  in  house  IDL  programs.  In 

165  each  tumor,  the  analysis  was  performed  on  data  from  3  slices,  each  2mm  thick. 

166 

167  The  MRSI  data  underwent  a  two  dimensional  spatial  Fourier  transform  and  the  spectral  grid  was 

168  superimposed  on  the  corresponding  T2-weighted  image  using  the  3DiCSI  software  package 

169  (courtesy  of  Truman  Brown,  Ph.D.,  Columbia  University).  The  voxels  within  the  tumor  were 
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170  identified  and  the  free  induction  decay  (FID)  from  each  tumor  voxel  was  extracted.  The  FID  was 

171  used  as  input  to  the  JMRUI  (Java  Magnetic  Resonance  User  Interface)  software  package  [28] 

172  and  the  lactate  resonance  was  fit  in  the  time  domain  using  AMARES  (Advanced  Method  for 

173  Accurate,  Robust  and  Efficient  Spectral  fitting).  The  2D-CSI  data  set  from  the  reference 

174  phantom  used  for  lactate  quantitation  (see  below)  was  processed  similarly.  For  quantitation, 

175  each  tumor  voxel  was  referenced  to  a  phantom  voxel  at  the  exact  same  location.  The  tumor  voxel 

176  lactate  peak  area,  the  reference  voxel  lactate  peak  area,  the  appropriate  T1  saturation  factors  and 

177  the  comparative  coil  loading  parameters  for  tumor  and  phantom  studies  were  used  to  quantify  the 

178  tumor  voxel  lactate  according  to  the  following  equation: 

179  Integration  area  of  unknown  X  Concentration  of  reference  standard  X  (  Vref/  Vtest) 

Unknown  concentration  =  -  - 

Integration  area  of  reference  standard  X  1-e(  ' 

180 

181  where  (Vref /Vtest)  is  ratio  of  the  voltage  difference  of  reference  sample  to  test  sample  and 

182  T1  is  the  spin-lattice  relaxation  of  lactate.  Lactate  in  the  doped  phantom  was  fully  relaxed.  In- 

183  vivo  T1  measurements  of  lactate  were  done  by  introducing  a  nonselective  180°  sine  pulse  to  the 

184  Sel-MQC  sequence  [29]  with  inversion  recovery  delays  from  10  ms  to  6  sec.  The  T1  value  for 

185  the  R3327-AT  tumor  (volume=1533  mm3)  was  1.26  ±  0.15  sec 

186 

187  To  evaluate  the  lactate  distribution  in  the  tumor  core  and  rim,  we  analyzed  lactate  in  both  the  rim 

188  (single  layer  of  outside  voxels  that  are  completely  in  the  tumor)  and  the  core  (all  other  voxels 

189  inside  the  tumor)  separately.  Those  voxels  were  also  used  to  define  rim  and  core  for  the  DCE- 

190  MRI  studies,  despite  the  difference  in  voxel  sizes  between  the  two  techniques. 

191 
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MR  Data  Statistics 

We  obtained  scans  for  each  rat  at  two  to  five  time  points.  Lactate  level  and  Akep  were 
determined  for  each  voxel  from  each  MRSI  image.  Lactate  was  recorded  as  a  numerical  value, 
or  as  zero  if  the  peak  was  not  detected.  Since  there  were  a  large  number  of  voxels  where  lactate 
was  not  detected,  especially  in  small  tumors,  either  treating  them  all  as  zero  or  excluding  them 
from  computations  would  introduce  bias.  Therefore,  we  assumed  that  the  log-lactate  levels  were 
normally  distributed,  and  that  any  voxel  without  detectable  lactate  was  considered  to  be  < 
1 .2mM  which  was  the  smallest  non-zero  value  observed  in  the  entire  experiment.  Lactate  and 
Akep  data  were  analyzed  using  the  generalized  linear  model. 

Measurements  were  obtained  for  varying  numbers  of  voxels  within  each  tumor  scan,  and  voxels 
were  classified  as  located  in  the  core  (C)  or  rim  (R)  of  the  tumor.  The  statistical  model  assumed 
that  1)  the  expected  value  of  lactate  followed  a  quadratic  function  of  tumor  volume,  2)  there  was 
random  variation  among  rats  in  a)  the  tumor  volume  at  which  lactate  reached  a  peak  and  b)  the 
magnitude  of  the  peak,  3)  there  was  a  proportional  difference  between  typical  values  for  core 
voxels  vs.  rim  voxels,  and  4)  the  variation  in  measurements  within  a  region  of  interest  followed  a 
censored  log-normal  distribution.  We  applied  a  non-linear  mixed  effects  statistical  model  using 
Proc  NLMixed  (SAS  version  9).  Two  rats  with  no  detectable  lactate  in  any  scan  yielded  no 
useful  information  and  were  excluded  from  the  analysis.  Analysis  of  Akep  values  was  more 
straightforward  because  there  were  very  few  zero  values.  Analysis  of  covariance  was  used  after 
applying  a  square  root  normalizing  transfonnation. 
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216 

217  HISTOLOGY 

218  Immunohistochemical  staining 

219  Four  rats  with  Dunning  R-3327-AT  tumors  were  sacrificed  following  serial  MR  examinations  for 

220  measurements  of  tumor  necrosis  and  hypoxia.  Since  initially  all  the  MR  experiments  were 

22 1  performed  serially  resulting  in  only  large  tumors  being  available  for  analysis,  an  additional  4  rats 

222  were  studied  by  MR,  each  at  a  single  time  point  and  sacrificed,  to  obtain  an  estimate  of  necrosis 

223  and  lactate  at  small  and  medium  tumor  volumes.  Six  additional  rats  were  sacrificed  at  various 

224  tumor  volumes  without  undergoing  MR  studies  specifically  to  provide  estimates  of  the  variation 

225  of  tumor  necrosis  with  volume 

226 

227  After  the  final  MRSI  experiments,  tumors  studied  by  histology  and  immunohistochemistry  were 

228  excised,  snap-frozen  in  isopentane/dry  ice,  embedded  in  cutting  medium  (OCT  4583,  Sakura 

229  Finetek)  and  8  pin-thick  contiguous  frozen  sections  from  multiple  locations  within  the  tumor 

230  obtained  using  a  cryostat  microtome  (Microm  International  GmbH).  Sections  were  fixed  in  4% 

231  paraformaldehyde  solution  for  10  minutes  at  room  temperature  then  blocked  in  Superblock-PBS 

232  (Pierce,  USA)  for  1  hour,  followed  by  either  immunofluorescence  staining  for  pimonidazole 

233  distribution,  or  haematoxylin/eosin  for  detennination  of  tumor  necrosis.  To  determine  the 

234  relationship  between  viable  tissue  volume  and  lactate,  the  necrotic  tumor  fraction  was  excluded 

235  for  these  calculations  by  estimating  the  percent  necrotic  area  on  5  histological  sections  taken 

236  from  separate  regions  of  each  tumor.  Pimonidazole  staining  was  carried  out  as  described 

237  previously  [30]  .  As  a  control  for  non-specific  antibody  binding,  some  sections  were  stained  with 

238  non-specific  mouse  IgG  conjugated  to  FITC. 


239 
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240  Images  of  pimonidazole  stained  tumor  sections  were  acquired  at  high  magnification  (xlOO)  using 

241  an  Olympus  BX40  fluorescence  microscope  (Olympus  America  Inc.,  Melville,  NY)  equipped 

242  with  a  motorized  stage  (Prior  Scientific  Instruments  Ltd.,  Cambridge,  UK)  and  an  appropriate 

243  filter.  The  individually  captured  image  frames  were  rendered  into  a  montage  of  the  entire  tumor 

244  section  using  ImagePro  software  (Media  Cybernetics  Inc,  Bethesda,  MD). 

245 

246  Histological  Data  Analysis 

247  For  determination  of  tumor  necrotic  fraction,  images  of  whole  tumor  sections  stained  with 

248  haematoxylin/eosin  were  examined  using  Photoshop  7.0  software.  Fractional  areas  of  necrosis 

249  were  determined  subjectively  for  each  section  by  two  investigators,  and  results  were  pooled. 

250  Sections  from  at  least  3  different  locations  within  a  single  tumor  were  used  to  determine  the 

25 1  mean  tumor  necrotic  fraction. 

252 

253  For  determination  of  mean  pimonidazole  binding  fraction,  immunofluorescence  images  of  anti- 

254  pimonidazole  images  were  converted  to  binary  form  using  a  threshold  determined  from  the 

255  staining  intensity  observed  in  tumor  sections  stained  with  the  control  antibody.  This  threshold 

256  remained  constant  for  all  tumors  analyzed.  The  percentage  of  tumor  area  with  staining  intensity 

257  greater  than  the  threshold  value  was  then  determined  using  Photoshop  7.0  software.  Sections 

258  from  at  least  3  different  locations  within  a  single  tumor  were  used  to  determine  the  mean  tumor 

259  pimonidazole-positive  fraction. 

260 
261 
262 

263 

264 

265 
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266  Experimental  Results 

267  In-vivo  lactate  measurement 

268  The  anaplastic,  faster-growing  R3327-AT  cell  line  had  an  in  vivo  tumor  doubling  time  (TDT)  of 

269  3.3  ±  0.17  days.  The  well-differentiated,  slower  growing  Dunning  H  had  a  doubling  time  of  22.8 

270  ±1.95  days  (Figure  1)  (p<0.001),  compatible  with  previous  studies  [20]. 

271 

272  Figure  2A  shows  a  series  of  T2  weighted  MR  images  of  5mm  thickness  co-registered  with  2D- 

273  CSI  lactate  spectra  obtained  on  the  Dunning  R3327  AT  tumor.  Lactate  peaks  in  the  tumor  region 

274  are  highlighted.  Lactate  was  not  detected  in  R3327-AT  tumors  studied  in  the  volume  range  110- 

275  265  mm  (group  I;  data  not  shown).  As  tumor  volumes  increased  to  the  range  of  478-649  mm 

276  (group  II),  the  average  voxel  lactate  concentration  detected  was  1.54  ±  0.54  mM  (Fig  3A).  In 

277  group  III  (tumor  volume  1035-1207  mm  ),  the  lactate  concentration  increased  to  4.79  ±  1.67 

278  mM.  In  Groups  IV  (1409-1708  mm3)  and  V  (1941-2453  mm3),  lactate  was  4.64  ±  1.26  mM  and 

279  2.71  ±  0.93  mM  respectively  (Fig  3A,  Table  1).  The  variation  of  average  lactate  concentration 

280  with  growth  is  shown  graphically  in  Figure  3 A  and  demonstrates  the  initial  increase  and 

28 1  subsequent  decrease  in  lactate  noted  with  tumor  growth.  Results  of  the  non-linear  mixed  effects 

282  analysis  showed  very  significant  curvature  in  the  relationship  between  lactate  and  tumor  volume 

283  (p=0.008)  namely,  an  initial  increase  followed  by  a  decrease  after  reaching  a  peak  that  varied 

284  among  rats  at  volumes  of  1048  +/-  467  (SD)  mm  .  There  was  a  significant  quadratic  or  peak 

285  effect  in  the  relationship  between  lactate  per  voxel  and  tumor  volume  (p<0.001)  in  the  R3327- 

286  AT  tumor  with  significant  variation  among  rats  in  the  volume  at  which  the  maximum  expected 

287  value  of  (log)  lactate  occurred. 

288 
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When  the  lactate  concentration  in  the  outer  voxels  of  the  tumor  (“rim”)  was  compared  to  all  the 
other  voxels  (“core”)  (Fig.  3B),  lactate  was  significantly  higher  in  the  core  compared  to  the  rim, 
independent  of  tumor  volume  (p<0.001).  The  lactate  concentration  vs  tumor  volume  curve  for 
the  core,  rim  and  whole  tumor  increases  to  a  maximum  followed  by  a  decrease  after  reaching  a 
peak  or  plateau  in  groups  III  and  IV.  This  was  characterized  by  a  quadratic  term  in  the  statistical 
model,  with  a  coefficent  that  was  significantly  less  than  zero,  reflecting  a  downturn  in  lactate 
concentration  with  volume  after  a  peak.  Inclusion  of  an  additional  parameter  in  the  statistical 
model  indicated  weak  evidence  (p=0.09)  that  on  average,  the  lactate  maximum  concentration 
occurred  at  larger  tumor  volumes  for  the  core  voxels  compared  to  the  rim. 

A  representative  set  of  2DCSI  spectra  of  Dunning  H  prostate  model  at  different  volumes  is 
represented  in  Figure  2B  (groups  I  to  V).  Lactate  was  not  detected  in  any  tumor,  regardless  of  the 
volume. 


DCE-MRI  analysis 
Dunning  R3327-AT 

The  average  Akep  for  group  I  tumors  was  1.59  ±  0.53  /min.  With  an  increase  in  tumor  size  (478 
-  649  mm  ),  Akep  decreased  to  0.65  ±  0. 12.  With  further  increases  in  tumor  size,  the  Akep  value 
for  group  III  and  IV  tumors  were  0.55  ±0.11  and  0.36  ±  0.04  respectively.  The  largest  tumors 
(group  V)  had  a  value  of  Akep  of  0.34  ±  0.05  (Fig.  4A;  Table  2A). 


We  analyzed  whether  the  rim  and  core  Akep  measurements  differed.  The  same  criterion  for 
classifying  the  rim  and  core  voxel  was  adopted  as  was  used  in  the  lactate  analysis  (Fig.  4B).  The 
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315  average  tumor  rim  and  core  Akep  value  of  the  lowest  tumor  volume  cohort  were  1 .6 1  ±  0.52 

316  /min  vs.  and  0.15  ±  0.14  /min  respectively,  but  only  one  of  the  4  animals  had  a  tumor  large 

317  enough  to  delineate  core  voxels  so  the  biological  significance  of  this  difference  is  unclear.  With 

318  the  increase  in  tumor  volume  to  478-649  mm  ,  the  tumor  rim  and  core  Akep  values  were  0.7 1  ± 

319  0.12  /min  vs.  0.27  ±  0.07  /min.  As  expected,  there  are  a  greater  number  of  voxels  in  the  rim 

320  than  in  the  core.  At  higher  tumor  volumes  (1035-1207  mm  and  1409-1708  mm  ),  the  rim  and 

321  core  Akep  values  decreased  to  0.63  ±  0.12  /min  vs  0.36  ±  0.12  /min  and  0.50  ±  0.07  /min  vs  0.15 

322  ±  0.03  /min  respectively.  At  the  largest  tumor  volumes  the  rim  and  core  Akep  values  were  0.47  ± 

323  .07/min  and  0. 14  ±  0.02/min  respectively.  The  rim  Akep  values  were  significantly  higher  than 

324  core  values  for  the  Dunning  R3327-AT  tumors  (p<0.001) 

325 

326  DH  tumors 

327  Slow  growing  Group  I  Dunning  H  tumor  had  an  average  Akep  value  of  1 .23  ±  0. 12/min  (Fig. 

328  5 A,  Table  2B).  This  decreased  to  1.06  ±  0.23/min,  0.84  ±  0.21/min  and  0.86  ±  0.08/min  in 

329  groups  II,  III  and  IV  tumors  respectively.  With  continued  tumor  growth,  the  whole  tumor  Akep 

330  values  decreased  (Group  V  =  0.46  ±  0.04/min).  We  subsequently  subdivided  the  analysis  to 

33 1  study  the  rim  and  core  as  we  had  done  previously  for  the  Dunning  R3327-AT  tumor  except  in 

332  the  Group  I  tumors  wherein  there  were  no  core  voxels  because  of  the  small  tumor  size.  The  rim 

333  and  core  Akep  values  both  decreased  with  growth  (p<0.001). 

334 

335  The  Dunning  H  exhibited  heterogenous  perfusion  as  manifested  by  significant  differences  in 

336  Akep  value  between  the  rim  and  core  (p<0.023).  The  Dunning  H  tumors  had  significantly  higher 

337  Akep  values  than  the  R3327-AT  tumors,  as  measured  in  the  whole  tumor  (p<0.001),  the  core 
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338  (p<0.001)  and  the  rim  (p<0.001).  Akep  declined  progressively  over  the  range  of  tumor  sizes  in 

339  the  Dunning  H  tumors  (p<0.001),  particularly  at  the  largest  tumor  volumes.  After  applying  an 

340  appropriate  normalizing  transformation  to  Akep  (e.g.,  cube  root),  the  rate  of  decrease  was  linear. 

341  This  decrease  was  also  noted  if  the  data  were  analyzed  separately  for  the  rim  and  core  (p<0.001). 

342 

343  Histology  and  Im munoh istoch em istry  Correlation 

344  The  Dunning  R3327-AT  tumors  have  necrosis  present  even  at  small  tumor  volumes  (Fig.  6a,  b, 

345  Fig.  7A).  In  contrast,  the  Dunning  H  tumors  had  no  observable  necrosis  even  at  volumes 

346  approaching  2300  mm  ,  (Figure  6c).  In  the  R33327-AT  tumor,  we  found  that  with  increasing 

347  tumor  volume,  the  necrotic  fraction  increased  significantly  (R  =  0.808  ;  p  <0.001)  (Fig,  7A).  A 

348  linear  correlation  between  the  concentration  of  lactate  measured  and  the  viable  tumor  volume 

349  was  noted  indicating  that  the  concentration  of  lactate  increases  with  tumor  growth  corrected  for 

350  viability  (R2  =  0.764;  p=0.005).  (Figure  7B) 

351 

352  Tumor  sections  sequential  to  those  used  for  histological  staining  were  assessed  for  the 

353  distribution  of  the  hypoxia  tracer  pimonidazole  by  immunofluorescence  staining  (Fig.  6(d),  (e), 

354  (f)).  We  observed  that  the  fractional  pimonidazole-positive  tumor  area  increased  with  tumor 

355  volume  (R“  =  0.720,  p  =  0.007).  Comparison  of  pimonidazole-positive  tumor  fraction  to  lactate 

356  values  also  showed  a  statistically  significant  positive  correlation  (R~  =  0.602,  p  =  0.023;  Fig. 

357  7C). 


16 


Lactate  imaging  and  DCE-MRI  in  prostate  tumors 

358  Discussion 

359  Tumor  lactate  kinetics  in  R332 7 AT  model 

360  Tumors  have  increased  lactate  levels,  due  to  increased  production  and/or  decreased  clearance.  Based  on 

361  this  study  of  two  prostate  tumor  models  with  different  growth  rates,  necrosis,  and  radiation 

362  sensitivity  [20],  we  hypothesize  that  lactate  is  a  potential  marker  of  prostate  tumor 

363  aggressiveness.  Lactate  concentrations  have  been  shown  to  have  prognostic  value  [13],  [14]  [31] 

364  [32]  in  other  tumors,  and  have  correlated  with  a  high  probability  of  metastasis  and  decreased 

365  patient  survival. 

366 

367  In  this  study,  lactate  was  not  detected  in  the  Dunning  R3327-AT  at  small  tumor  volumes  (1 10  - 

368  265  mm  ).  At  the  lowest  tumor  volumes,  the  tumor  appears  to  be  effectively  and  uniformly  well 

369  perfused,  based  on  the  measurements  of  both  lactate  and  Akep,  suggesting  adequate  oxygen  and 

370  nutrient  delivery  to  meet  the  metabolic  requirements  of  the  tumor.  As  the  R3327-AT  tumor 

371  volume  increased,  lactate  was  detected  and  was  associated  with  a  concomitant  decrease  in  Akep 

372  and  metabolic  and  perfusion  heterogeneity,  suggesting  that  the  tumor  has  evolved  into  a  more 

373  metabolically  deprived  state  associated  with  poorer  perfusion.  The  perfusion  and  metabolic 

374  deficits  become  more  pronounced  as  tumor  growth  continues  and  is  compatible  with  induction  of 

375  hypoxia  (Fig.  6e)  and  the  tumor  being  forced  to  undergo  increasing  anaerobic  glycolysis  with 

376  increased  lactate  production. 

377 

378  In  contrast,  in  the  Dunning  H  lactate  is  not  detected  and  there  is  a  significant  drop  in  Akep 

379  primarily  at  the  highest  tumor  volume  range.  The  Akep  values  were  approximately  constant 

380  between  1000  and  1700  mm  .  These  relatively  constant  values  of  Akep  in  this  volume  range  are 
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381  noted  both  in  the  rim  and  the  core  of  the  tumor  (Fig.  5,  Table  2B).  These  data  explain  the  lack 

382  of  necrosis,  hypoxia  and  lactate  in  this  model.  We  hypothesize  that  the  relative  absence  of  lactate 

383  in  the  Dunning  H  is  due  to  the  better  perfusion,  minimizing  hypoxia,  anaerobic  glycolysis,  and 

384  necrosis.  The  absence  of  lactate  at  very  large  tumor  volumes  wherein  perfusion  has  decreased 

385  (although  the  decrease  in  Akep  is  not  as  large  as  in  the  Dunning  R3327-AT)  may  be  due  to  the 

386  slower  growth  rate  of  the  Dunning  H  which  would  place  a  lower  demand  on  energy  requirements 

387  and  therefore  lower  glucose  requirements.  A  major  cause  of  the  differences  in  the  tumor  models 

388  may  be  ascribed  to  differences  in  perfusion  as  measured  by  Akep  which  is  significantly  higher  in 

389  the  Dunning  H  than  in  the  R3327-AT. 

390 

391  Lactate  was  detected  prior  to  the  presence  of  significant  histologic  necrosis  (Fig.  7A)  in  the 

392  R3327-AT  tumor  indicating  that  the  lactate  is  not  solely  due  to  the  presence  of  necrotic  tissue. 

393  Figures  3A  and  7A  demonstrate  that  lactate  decreases  at  larger  tumor  volumes  where  necrosis 

394  becomes  more  significant,  indicating  that  lactate  is  produced  in  viable  tissue.  In  addition,  there 

395  is  a  linear  increase  in  lactate  as  viable  tumor  volume  increases,  further  indicating  that  lactate 

396  occurs  independently  of  necrosis. 

397 

398  Previous  in  vitro  studies  have  indicated  that  lactate  is  elevated  in  prostate  cancer  [33],  [34]  [35], 

399  Using  hyperpolarized  13C  spectroscopic  imaging,  Chen  et  al  [36]  studied  the  transgenic 

400  adenocarcinoma  of  mouse  prostate  (TRAMP)  and  found  that  after  infusion  of  13C  labeled 

401  pyruvate,  lactate  could  be  detected  in  both  primary  and  metastatic  prostate  tumors.  Albers  et  al 

402  [37]  correlated  lactate  levels  with  a  histologic  index  developed  by  the  authors  and  suggested  that 

403  lactate  levels  increased  as  the  tumor  grade  increased,  as  measured  by  their  histologic  index. 
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404 

405  Numerous  recent  studies  have  shown  that  tumor  hypoxia  and  altered  metabolism  are 

406  interconnected.  Hypoxia  is  known  to  enhance  glycolytic  activity,  and  likely  works  via  the  HIF1- 

407  a  pathway.  HIFl-a  is  increased  under  hypoxic  conditions  and  increases  the  expression  of 

408  numerous  genes  related  to  glucose  metabolism  including  glucose  transporters  [38],  glycolytic 

409  enzymes  [39]  and  also  suppresses  the  activity  of  the  citric  acid  cycle.  Recent  studies  have 

410  demonstrated  that  HIF1-  a  causes  overexpression  of  pyruvate  dehydrogenase  kinase- 1  (PDK-1) 

411  [40-41]  which  lowers  pyruvate  dehydrogenase  (PDH)  activity.  Loss  of  PDH  activity  suppresses 

412  conversion  of  pyruvate  to  acetate  thereby  diminishing  citric  acid  cycle  activity,  oxidative 

413  phosphorylation,  oxygen  consumption  and  formation  of  reactive  oxygen  species  (ROS).  In  vivo 

414  lactate  measurements,  could  be  used  to  monitor  the  activity  of  targeted  drugs  aimed  at  inhibiting 

415  various  metabolic  pathways,  and  used  to  determine  an  appropriate  dose  and  validating  it  in 

416  patients  by  monitoring  lactate  and  changes  in  its  concentration. 

417 

418  In  addition  to  verifying  that  a  drug  is  active  at  the  target  site,  the  in  vivo  measurement  of  lactate 

419  in  prostate  cancer  could  yield  both  a  prognostic  marker  as  has  been  shown  previously  [37]  and 

420  would  potentially  be  useful  for  “dose  painting”  by  detecting  areas  of  potential  hypoxia  and/or 

421  radiation  resistance  [42-43].  One  of  the  major  advantages  to  obtaining  localized  lactate  data 

422  instead  of  measurements  from  the  whole  tumor  is  the  ability  to  delineate  the  heterogenous  state 

423  of  the  tumor  as  shown  in  this  study  and  potentially  applying  it  for  local  therapies. 

424 

425  Prostate  cancer  is  a  very  important  tumor  to  study  because  of  its  high  incidence  associated  with  a 

426  high  survival  rate,  and  often  indolent  course.  Studies  have  demonstrated  that  many  patients  with 
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427  prostate  cancer  do  not  need  treatment  since  it  does  not  alter  the  outcome  of  the  disease  [44]. 

428  There  are  currently  no  reliable  methods  of  selecting  these  low  risk  patients  and  this  is  an  area  of 

429  intense  ongoing  research  [45].  The  current  study  suggests  that  measurements  of  lactate  in 

430  prostate  tumors  might  be  an  additional  appropriate  biomarker  for  investigations  and  inclusion 

43 1  with  other  clinical  parameters  into  nomograms.  Since  proton  MRSI  data  of  the  prostate  is  often 

432  obtained  in  the  course  of  staging  workup  after  diagnosis,  these  data  are  readily  translated  to  the 

433  clinic.  Lactate  editing  methodologies  do  not  require  contrast  agents  or  isotopic  injection  and  are 

434  likely  to  be  sensitive  enough  to  detect  lactate  in  patients  [46-48].  The  use  of  13C  hyperpolarized 

435  MRSI  [37]  would  significantly  enhance  the  technique,  if  it  receives  clinical  approval. 

436 

437  Estimation  of  tumor  hypoxic  percentage  (using  pimonidazole  immunofluorescence  staining) 

438  showed  a  positive  correlation  with  increasing  lactate  in  R3327-AT  tumors  (Fig.  7C)).  In  keeping 

439  with  previous  studies  [2 1],  we  were  unable  to  detect  either  significant  pimonidazole  uptake  or 

440  lactate  in  any  of  the  Dunning-H  tumors  examined.  Taken  together,  these  observations  are 

441  supportive  of  the  hypothesis  that  hypoxic  tumor  tissue  is  a  significant  source  of  tumor  lactate 

442  production.  However,  as  the  pimonidazole-positive  tumor  area  also  increases  with  tumor  volume 

443  in  the  R3327-AT,  it  is  also  possible  that  increasing  tumor  volume  with  inadequate  angiogenesis 

444  is  a  contributing  cause  to  lactate  production. 

445 

446 

447 

448 


449 
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450  Conclusion 

45 1  Lactate  is  not  present  in  the  slowly  growing  Dunning  H  but  is  detected  in  a  heterogenous  pattern 

452  in  the  fast  growing  R3327-AT  tumor,  suggesting  that  in  vivo  lactate  measurements  could  be  a 

453  marker  of  aggressive  tumors.  Based  on  the  tumor  heterogeneity  in  lactate  concentration,  we 

454  hypothesize  that  localized  areas  of  aggressive,  hypoxic  tumors  could  be  potentially  delineated  for 

455  localized  therapy.  The  core  of  the  tumor  is  prone  to  accumulate  more  lactate  than  the  rim  of  the 

456  tumor,  probably  since  the  core  of  the  tumor  would  be  expected  to  have  a  less  adequate  oxygen 

457  and  blood  supply.  This  observation  is  supported  by  the  DCE-MRI  study  where  the  tumor  core  is 

458  significantly  less  well  perfused  than  its  rim  in  the  R3327-AT  tumor. 

459 
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468  Figure  Captions 

469  Figure  1.  Tumor  volume  (mean  ±  SEM)  measurements  of  Dunning  R3327  AT  (A)  and  Dunning 

470  H  (B)  tumors.  Note  different  time  scales 

471 

472  Figure  2  A.  2DCSI  lactate  spectra  of  R3327  AT  prostate  tumor  model  as  a  function  of  tumor 

473  volume,  1:613  mm3,  II:  1054  mm3.  III:  1455  mm3  and  IV:  2207  mm3.  T2  weighted  saggital  MR 

474  tumor  image  is  co-registered  with  2DCSI  lactate  spectra.  The  peak  intensity  of  the  lactate 

475  resonance  is  heterogenous,  even  within  a  specific  tumor. 

476 

477  Figure  2B.  2DCSI  lactate  spectra  of  Dunning  H  prostate  tumor  model  as  a  function  of  tumor 

478  volume,  A:  246  mm3  B:512  mm3,  C:1075  mm3,  D:1382  mm3  and  E:  2255  mm3.  T2  weighted 

479  saggital  MR  tumor  image  is  co-registered  with  2DCSI  lactate  spectra.  Note  the  absence  of 

480  lactate  in  any  of  the  voxels 

481 

482  Figure  3A.  Measurement  of  average  lactate  in  R3327  AT  tumor  model  as  a  function  of  tumor 

483  volume.  Tumors  in  the  volume  range  of  1 10-  265  mm3  did  not  show  lactate 

484 

485  Figure  3B.  Evaluation  of  tumor  core/rim  lactate  of  R3327  AT  tumor  model  as  function  of  tumor 

486  model  (solid  filled  bar  plots=  core  lactate;  empty  bar  plots=  rim  lactate).  Lactate  initially 

487  increases,  shows  a  plateau,  and  then  decreases  in  the  whole  tumor,  rim,  and  core 

488 

489  Figure  4A.  Whole  tumor  average  Akep  (/min)  of  Dunning  R3327  AT  tumor  model  as  a  function 

490  of  tumor  volume.  Akep  decreases  with  tumor  volume 


22 


Lactate  imaging  and  DCE-MRI  in  prostate  tumors 

491  Figure  4B.  Evaluation  of  tumor  core/rim  Akep  (/min)  of  Dunning  R3327  AT  tumor  model  at 

492  different  tumor  volumes. 

493 

494  Figure  5  A.  Whole  tumor  average  Akep  (/min)  of  Dunning  H  tumor  as  a  function  of  tumor 

495  volume. 

496 

497  Figure  5B.  Evaluation  of  tumor  core/rim  Akep  (/min)  of  Dunning  H  tumor  model  over  the 

498  volume  range  studied. 

499 

500  Figure  6.  H  &  E  staining  of  (a)  Dunning  R3327  AT  (118  mm3),  (b)  Dunning  R  3327-AT  (2207 

501  mm  )  and  (c)  Dunning  H  (  1949  mm  ).  Corresponding  slides  showing  pimonidazole  staining  are 

502  displayed  in  (d),  (e)  and  (f). 

503 

504  Figure  7A.  A  graph  of  percentage  of  necrosis  as  a  function  of  tumor  volume  demonstrates  a 

505  significant  linear  correlation  between  these  two  parameters. 

506 

507  Figure  7B.  A  graph  of  the  variation  of  tumor  lactate  vs.  viable  tumor  volume  demonstrating  a 

508  significant  linear  correlation.  Lactate  is  noted  to  increase  with  tumor  growth,  when  corrected  by 

509  eliminating  the  necrotic  tissue 

510 

511  Figure  1C.  Graph  of  pimonidazole  positive  fraction  as  a  function  of  measured  lactate  showing  a 

512  linear  correlation  between  pimonidazole  and  lactate  concentration.  The  data  suggest  that  lactate 

513  is  a  potential  surrogate  for  measuring  tumor  hypoxia  in  this  model. 
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5 14  Table  1  Average  lactate  concentration  and  numbers  of  CSI  voxels  in  the  Dunning  R-3327-AT 

515  tumor 

516 


Whole  tumor 

Core 

Rim 

Tumor 

Group 

Tumor  volume 
(mm3) 

Lactate 

(mM/L) 

Avg# 

voxels 

Lactate 

(mM/L) 

Avg# 

voxels 

Lactate 

(mM/L) 

Avg# 

voxels 

I 

110-265 

0 

7 

0 

1 

0 

6 

II 

478-649 

1.54 

13 

1.79 

3 

1.46 

10 

III 

1035-1207 

4.79 

17 

6.95 

6 

3.68 

11 

IV 

1409-1708 

4.64 

20 

7.34 

8 

3.02 

12 

V 

1941-2453 

2.71 

24 

4.78 

10 

1.40 

14 
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519  Table  2A  Average  Akep  values  of  Dunning  R  3327  AT  tumor 

520 


Tumor 

volume 

(mm3) 

Tumor 

group 

Whole 

tumor 

Akep 

Rim 

Akep 

Average 
No.  of 
rim  CSI 
voxels 

Core 

Akep 

Average 
no.  of 
core  CSI 
voxels 

110-265 

I 

1.59±0.53 

1.61±0.52 

6 

0.15±0.14 

1 

478-649 

II 

0.65±0.12 

0.71±0.12 

10 

0.27±0.07 

3 

1035-1207 

III 

0.55±0.1 1 

0.63±0.12 

11 

0.36±0.12 

6 

1409-1708 

IV 

0.36±0.04 

0.50±0.07 

12 

0.15±0.03 

8 

1941-2453 

V 

0.34±0.05 

0.47±0.07 

14 

0.14±0.02 

10 
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52 1  Table  2B  Average  Akep  values  of  Dunning  H  tumors 

522 


521 


Tumor 

volume 

(mm3) 

Tumor 

group 

Whole 

tumor 

Akep 

Rim 

Akep 

Average 
no.  of 
Rim  CSI 
voxels 

Core 

Akep 

Average 
no.  of 
Core 
CSI 
voxels 

227-259 

I 

1.23±0.12 

1.23±0.12 

7 

0 

430-546 

II 

1.06±0.23 

1.07±0.23 

9 

0.76±0.27 

1 

1000-1236 

III 

0.84±0.21 

0.85±0.21 

11 

0.77±0.23 

3 

1382-1715 

IV 

0.86±.08 

0.90±0.04 

13 

0.75±0.19 

7 

1886-2255 

V 

0.46±.04 

0.53±0.05 

12 

0.38±.06 

11 
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Figure  2 
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Figure  5A  Dunning  H 
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529  Figure  6 
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